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DETECTION  PERFORMANCE  SIMULATION 
OF  A  DICKE-TYPE  CFAR  IN 
DISTRIBUTED  CLUTTER 


J.  N.  Bucks sm 


ABSTRACT 


^A  Dicke-type  CFAR  processor  is  one  which  hard  limits  the  radar 
IF  or  bipolar  video  signal  prior  to  pulse  compression.  The  constant  false 
alarm  rate  (CFAR)  properties  of  this  processor  in  white  receiver  noise 
are  well  known.  When  the  background  against  which  targets  must  be 
detected  is  due  to  clutter  returns,  this  CFAR  property  is  destroyed. 


This  report  summarises  results  of  computer  simulation  of  a 

Dicke-type  CFAR  in  a  background  consisting  of  distributed  clutter  returns. 

Clutter  model  parameters  are  chosen  to  approximate  the  statistics  of 

sea,  rain,  and  land  clutter.  The  simulations  show  the  increase  in  pro- 

n 

bability  of  false  alarm  (P  ')  in  these  backgrounds  when  the  threshold 

FA  y_. 

settings  are  chosen  to  provide  the  desired  P^  in  receiver  noise.  In 
addition,  the  probability  of  target  detection  as  a  function  of  signal- 
to-background  ratio  is  determined.  Values  are  determined  for  the 
coefficients  in  a  parametric  representation  of  these  results. 


A  detailed  documentation  of  the  simulation  software  Is  included 


Is  incl 


as  an  Appendix. 
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1.0 


INTRODPCTION 


A  Dicke-type  CFAR  processor  is  one  which  herd  limits  the  radar 
IF  or  bipolar  video  signal  prior  to  pulse  compression.  The  constant  false 
alarm  rate  (CFAR)  properties  of  this  processor  in  white  receiver  noise 
are  well  known  [1].  When  the  background  against  which  targets  oust  be 
detected  is  due  to  clutter  returns,  this  CFAR  property  is  destroyed. 

This  report  summarizes  results  of  computer  simulation  of  a 
Dicke-type  CFAR  in  a  background  consisting  of  distributed  clutter  returns. 
Clutter  model  parameters  are  chosen  to  approximate  the  statistics  of 
sea,  rain,  and  land  clutter.  The  simulations  show  the  increase  in  pro¬ 
bability  of  false  alarm  (Py^)  in  these  backgrounds  when  the  threshold 

settings  are  chosen  to  provide  the  desired  P  in  receiver  noise.  In 

fA 

addition,  the  probability  of  target  detection  as  a  function  of  signal-to- 
background  ratio  is  determined.  Values  are  determined  for  the  coefficients 
in  a  parametric  representation  of  these  results. 


2.0  PROCESSOR  AMD  BACKGROUND  DESCRIPTIONS 

2.1  Waveform  and  Signal  Processor 

Figure  1  depicts  the  receive  signal  processing. 

The  transmitted  waveform  consists  of  a  phase-coded  pulse-  com¬ 
pression  waveform.  It  Is  formed  by  a  sequence  of  contiguously  transmitted, 
equal  amplitude  pulses  (or  "chips),  each  of  which  has  a  phase  of  either 
04  or  180* ,  determined  by  a  maximal-length  pseudo-random  sequence  generator. 
This  waveform  Is  transmitted  sequentially  on  one  or  more  frequency  di¬ 
versity  channels. 

On  receive,  the  IF  signal  (on  a  particular  frequency  diversity 
channel)  Is  first  converted  to  I  and  Q  bipolar  video  by  filters  matched 
to  the  individual  chips.  The  bipolar  video  is  hard-limited  to  +/-  one, 
and  sampled  at  the  chip  rate.  The  hard-limited  samples  are  then  pulse 
compressed  by  cross-correlation  against  the  stored  pseudo-random  sequence 
of  +/-  ones.  The  pulse-compressed  I  and  Q  are  then  square- law  envelope 
detected,  and  added  to  the  corresponding  quantities  from  other  diversity 
channels.  This  sum  forms  the  test  statistic  for  threshold  comparison. 


2.2  Background  Models 

Four  types  of  background  are  of  interest: 

1.  Receiver  noise 

2.  Distributed  Sea  Clutter 

3.  Distributed  Rain  Clutter 

4.  Distributed  Land  Clutter 

Receiver  noise  is  easily  (and  accurately)  described  as  Gaussian 
noise,  which  is  spectrally  flat  (white)  over  the  receiver  bandwidth.  Thus, 
samples  of  receiver  noise  are  independent  when  separated  in  time  by  roughly 
the  inverse  of  the  receiver  bandwidth. 
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•Limiting  Pulse  Compression  Receiver  with  Incoherent 
i-to-Pulse  Integraion 


The  radar  clutter  response  requires  a  relatively  more  compli¬ 
cated  model.  The  clutter  return  from  a  given  resolution  cell  is  assumed 
to  be  the  sum  of  contributions  from  a  large  number  of  independent  "point" 
clutter  scatterers  within  the  cell.  When  the  number  of  such  scatterers 
is  large,  and  when  the  sum  is  not  dominated  by  a  small  subset  of  scatterers, 
samples  of  the  radar  clutter  signal  produced  by  the  Isolated  resolution 
cell  on  successive  pulses  may  be  modeled  as  a  Rayleigh  random  process, 
whose  power  is  proportional  to  the  number  and  strength  of  scatterers 

within  the  cell.  On  any  given  cell,  however,  the  mean  power  of  this 
Rayleigh  process  is  apt  to  differ  quite  significantly  from  other  cells 
within  the  local  vicinity.  This  cell  to  cell  variation  of  the  local 
average  cross  section  is  more  accurately  described  by  the  log-normal  or 
Weibul  families  of  distributions. 

This  leads  to  the  following  algorithm  for  simulating  the 
clutter  return  from  a  group  of  range  cells  on  a  sequence  of  radar  pulses. 

A  random  number  Y  is  drawn  for  each  cell  from  a  log-normal  or  Weibul 
distribution  to  represent  the  average  clutter  cross-section  occurring 
in  that  cell.  For  each  cell,  a  sequence  of  clutter  amplitudes  is  gen¬ 
erated  by  sampling  a  narrowband  complex  Gaussian  process  whose  mean 
power  is  the  average  cross-section  from  that  cell,  and  whose  sample-to- 
sample  correlation  corresponds  to  the  desired  pulse-to-pulse  fluctuation 
spectrum  of  the  clutter. 

Mathematically,  the  complex  clutter  amplitude  is  obtained  as 
the  product  of  two  random  variables: 


where 


and 


/7J  X 


ij 


is  the  complex  (in-phase  and  quadrature)  clutter  amplitude 
on  the  i**4  range  bin,  pulse, 

Y^  is  the  mean  clutter  cross-section  on  the  1th  range  cell, 
which  follows  a  log-normal  or  Welbull  distribution. 


is  a  unit  mean  power,  complex  Gaussian  random  variable. 
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The  complex  clutter  onplltude  is  a  model  for  the  complex  bipolar  video 
clutter  responae  of  a  linear,  simple  pulse  (no  pulse  compression)  radar 
with  pulse  length  equal  to  the  chip  length  of  the  pseudo-noise  code. 
That  is,  in  Figure  1,  represents  the  I,  Q  samples  just  prior  to  the 
hard  limiters,  when  a  unit  compression  ratio  signal  ft^l)  is  utilized: 

I  -  real  (C^) 

Q  -  lmag  (C^) 


When  a  pulse  compression  waveform  is  used, 
by  convolution  in  the  range  dimension  of  C 


the  I,  Q  samples  are  obtained 
^  with  the  transmit  waveform  samples 


Table  1  summarizes  the  specific  distributions  used  for  simulating  the 
cell-to-cell  variations  (Y^)  in  mean  clutter  cross-section,  for  each 
of  three  types  of  clutter.  For  sea  and  rain  clutter,  is  drawn  from 
a  log-normal  distribution.  That  is. 


Y 


e 


z 


where  z  is  normally  distributed.  The  distribution  for  Y  can  be  written 


where 

Y^q  is  the  median  value  of  Y 
and  S  is  the  standard  deviation  of  ln(Y/Y^). 

The  average  (Y)  and  median  (Y^q)  values  of  Y  are  related  by 

Y  -  Y50  exp{S2/2> 
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For  purposes  of  simulation,  the  average  value  (Y)  of  the  mean  clutter 
cross-section  (Y)  is  constrained  to  unity,  i.e.,  unit  average  mean  clutter 
cross-section.  The  remaining  parameter  in  the  distribution  is  the  log- 
standard  deviation  S.  Expressed  in  dB  units,  this  is 

SDB  -  /  E[10  logl0(Y/Y50)]2 
-  4.34  S 

For  sea  clutter,  SDB  is  6  dB,  and  for  rain  clutter,  3  dB.  The  mean  to 
median  ratios  are  thus 

7/t5o  -  «p  ji/2  <f?|/J 

12.60  ,  sea  clutter 

1.27  ,  rain  clutter 


For  land  clutter,  is  drawn  from  a  Welbull  distribution.  A 
unit  mean  tfeibul  probability  density  is 


c  c 

PyCY)  -  c  Yc*1tr(i  +  £)]  expi-rmi  +  £)]  ) 


The  mean-to-median  ratio  for  this  distribution  is 


Vh0- 


r«  + 


(ln2) 


1/C 


The  value  of  C  is  0.3,  for  which 


3.0 


SIMULATION  RESULTS 


Detection  performance  of  the  processor  described  in  section 

2.1  is  simulated  in  this  section.  The  pulse  compression  ratio  in  all 
cases  is  127:1.  Results  are  presented  for  1,  2,  and  4  diversity  channels. 
A  non-fluctuating  (Sverling  0)  target  is  assumed. 

The  simulations  entail  three  basic  steps.  First,  the  detec¬ 
tion  thresholds  are  selected  to  provide  the  desired  in  receiver 
noise.  Using  these  thresholds,  the  actual  P^'s  are  then  obtained  in 
each  of  the  three  types  of  clutter.  Finally ,  the  detection  probability 
is  determined  as  a  function  of  si gnal-to-back ground  ratio,  in  each 
clutter  type.  This  functional  dependence  is  then  summarized  by  a  four 
parameter  smooth  curve  which  is  fit  to  the  simulation  results  using 
the  technique  discussed  in  [2]. 

3.1  Threshold  Selections  in  Receiver  Noise 

Figure  2  summarizes  the  false  alarm  behavior  of  the  processor 
in  white,  Gaussian  receiver  noise,  as  determined  by  Monte  Carlo  simu¬ 
lation.  The  figure  is  based  on  103  trials.  It  was  found  empirically 
that  a  nearly  linear  relationship  existed  between  the  threshold  level 
in  dB  (10  log..  T,  where  T  is  the  threshold  value  in  Figure  1)  and 

■*-u  _4 

l°gl0  ^“log10  ^PFA^’  OVer  a  range  of  Pfa  from  -63  to  *10  *  This 
law  was  assumed  to  continue  for  P^  below  10  ,  so  that  the  threshold 

settings  in  Table  2  could  be  extrapolated  from  the  simulation  data. 

3.2  Probability  of  False  Alarm  in  Non-Gaussian  Clutter 

When  the  threshold  settings  for  receiver  noise  are  used  in 
other  background  such  as  distributed  clutter,  the  probabilities  of  false 
alarm  increase  by  orders  of  magnitude.  Figures  3  to  5  summarize  the 
false  alarm  performance  of  the  processor  for  the  distributed  sea,  rain, 
and  land  clutter  models  of  section  2.2.  The  figures  are  based  on  Monte 
Carlo  simulation  with  103  trials. 
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Table  2  -  Threshold  Settings  In  Receiver  Noise 


N  -  1 

N  -  2 

N  -  4 

32.1  dB* 

33.3 

35.0 

33.3 

34.2 

35.7 

34.2 

34.9 

36.2 

35.0 

35.4 

36.7 

35.6 

35.9 

37.0 

36.2 

36.3 

37.4 

♦tabulated  quantity  Is  10  log1Q  T,  where  T  is  the  threshold 
setting. 


3.3  Probability  of  Detection.  Non-Fluctuating  Target 

Detection  performance  for  each  background  types  was  simulated 
using  a  non-fluctuating  (Swerling  0)  target  model,  and  200  trials.  The 
threshold  settings  from  Table  2  were  used.  The  results  are  plotted 
in  Figures  6  to  17,  as  curves  of  probability  of  detection  (Pg)  vs. 
signal-to-back ground  ratio  (SBR).  When  the  background  is  receiver  noise, 
signal-to-back ground  ratio  is  defined  as  signal  energy  to  noise  spec¬ 
tral  density  ratio.  When  the  background  is  distributed  clutter,  signal- 
to-background  is  defined  as  the  ratio  of  average  target  cross-section 
to  average  mean  clutter  cross-section  per  range  bin.  In  each  figure, 
the  curves  are  keyed  to  the  false  alarm  probability  obtained  for  each 
threshold  in  the  particular  background  type. 

The  smooth  curves  connecting  the  simulation  data  points  were 
obtained  by  application  of  the  parametric  curve  fitting  techniques  of  [2]. 
The  coefficients  for  all  cases  in  Figures  6  to  17  are  summarized  in 
Table  3. 
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Figure  11 

PERFORMANCE  IN  SEA  CLUTTER;  M“t27;  N«2 
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SIGNAL  TO  BACKGROUND  COB) 

Figure  12 

PERFORMANCE  IN  RAIN  CLUTTER,  M“127;  N-2 


SIGNAL  TO  BACKGROUND  COB) 

Figure  13 

PERFORMANCE  IN  LAND  CLUTTER;  M-127;  N-2 
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Figure  14 
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Table  3  -  (continued) 
(b)  Sea  Clutter 


2.63 

X 

10'2 

9.86 

x  IO”3 

1.51 

25.85 

1.51 

X 

io"2 

1.69 

x  IO*2 

1.74 

10.8 

9.8 

X 

io“3 

3.45 

x  10~2 

2.14 

2.95 

6.5 

X 

io'3 

3.11 

x  10"2 

2.58 

2.72 

4.5 

X 

10~3 

2.51 

x  10“2 

2.70 

3.26 

3.0 

X 

io“3 

2.74 

x  10“2 

2.94 

1.96 

4.17 

X 

10~2 

8.15 

x  IO-2 

1.86 

3.71 

2.88 

X 

io'2 

8.60 

x  10“2 

2.12 

2.43 

2.04 

X 

10~2 

8.36 

x  10~2 

2.47 

2.13 

1.62 

X 

10~2 

6.43 

x  IO”2 

2.39 

2.42 

1.27 

X 

io'2 

6.86 

x  IO*2 

2.73 

1.82 

9.8 

X 

io'3 

6.71 

x  IO*2 

3.07 

1.59 

6.46 

X 

io"2 

1.51 

x  10'2 

2.30 

492 

5.01 

X 

io"2 

1.52 

x  IO'1 

2.94 

2.48 

4.07 

X 

IO'2 

1.23 

x  IO'1 

2.95 

2.93 

3.63 

X 

10‘2 

1.09 

x  IO'1 

3.03 

3.03 

2.95 

X 

io'2 

1.17 

x  IO'1 

3.29 

2.06 

2.63 

X 

io"2 

9.23 

x  IO'2 

3.15 

2.71 

II 


N 

1 


2 


2 


Table  3  -  (continued) 
(c)  Rain  Clutter 


I 

1.14 

FA 

x  10"2 

4.97 

A 

X 

io"3 

B 

1.73 

C 

136 

2.7 

x  IO-3 

1.50 

X 

io"2 

2.10 

22.3 

8.0 

x  10"3 

2.74 

X 

io'2 

2.54 

5.60 

2.2 

x  10"3 

2.87 

X 

io'2 

2.66 

3.26 

7.3 

x  10-4 

2.94 

X 

IO'2 

3.06 

2.45 

2.1 

-A 

x  10 

2.79 

X 

IO'2 

3.33 

1.95 

2.34 

x  10“2 

5.35 

X 

io-2 

1.59 

5.46 

9,8 

x  10‘3 

6.17 

X 

IO-2 

1.96 

4.14 

3.8 

x  10“3 

5.43 

X 

io-2 

2.29 

4.64 

1.9 

x  10"3 

5.0 

X 

io'2 

2.46 

4.50 

9.0 

x  IO-4 

6.01 

X 

IO"2 

2.84 

2.58 

3.9 

x  10"4 

5.57 

X 

io-2 

2.99 

2.44 

4.47 

x  10"2 

1.02 

X 

io-1 

2.24 

6.69 

2.29 

x  10"2 

9.10 

X 

io’2 

2.53 

7.01 

1.32 

x  10"2 

1.18 

X 

io-1 

2.99 

3.33 

9.10 

x  10"3 

1.30 

X 

10_1 

3.23 

2.13 

4.70 

x  10'3 

1.31 

X 

IO’1 

3.41 

1.55 

3.00 

x  10"3 

9.94 

X 

io-2 

3.30 

2.30 
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Table  3  -  (continued) 
(d)  Land  Clutter 


N 

PFA 

A 

B 

C 

1 

2.82 

x  IO-2 

1.71 

X 

IO-1 

2.55 

.596 

2.00 

x  10-2 

1.16 

X 

10_1 

2.52 

.688 

1.58 

x  IO-2 

8.29 

X 

10'2 

2.62 

.749 

1.26 

x  10-2 

6.88 

X 

IQ*2 

2.91 

.683 

1.05 

x  10‘2 

6.44 

X 

IQ'2 

3.20 

.598 

8.80 

x  10“3 

6.31 

X 

10-2 

4.17 

.410 

2 

4.47 

x  10“2 

5.95 

X 

10~2 

1.97 

6.46 

3.39 

x  10"2 

7.75 

X 

10”2 

2.12 

3.03 

2.69 

x  10‘2 

8.48 

X 

10‘2 

2.36 

2.05 

2.29 

x  10~2 

1.05 

X 

10'1 

2.75 

1.17 

1.96 

x  10~2 

9.48 

X 

10-2 

2.91 

1.08 

1.67 

x  10“2 

8.6 

X 

io-2 

2.91 

1.13 

4 

6.46 

x  10”2 

2.94 

X 

10_1 

2.65 

1.11 

5.37 

x  10"2 

2.25 

X 

io"1 

2.69 

1.18 

4.57 

x  IO-2 

2.24 

X 

io-1 

3.22 

9.60 

4.17 

x  10-2 

1.85 

X 

IO'1 

3.11 

1.10 

3.63 

x  10-2 

1.65 

X 

IO*1 

3.17 

1.15 

3.31 

x  10-2 

1.51 

X 

10"* 

3.26 

1.11 

4.0  SIMULATION  TECHNIQUES 

The  simulations  were  implemented  in  FORTRAN  IV  on  a  Prime  350 
computer.  This  section  summarizes  the  functional  organization  of  the 
simulation  program,  and  discusses  a  few  simplifications  utilized  in  the 
modeling.  A  more  detailed  description  of  the  program  modules  is  included 
as  Appendix  A. 

4.1  Functional  Description  of  Simulation  Program 

Figure  18  shows  the  functional  organization  of  the  simulation 
program,  at  the  top  level  of  module  hierarchy.  The  first  major  operation 
of  the  simulation  is  the  generation  of  a  numerical  representation  of  the 
transmit  waveform.  This  is  done  with  a  call  to  the  subroutine  PNCODE, 
which  returns  a  sequence  of  plus  and  minus  ones  generated  from  a  maximal 
length  shift  register.  This  sequence  represents  the  sequence  of  amplitudes 
(or  equivalently  phases)  transmitted  on  each  chip  in  the  phase  coded 
waveform. 


The  second  step  is  the  generation  of  the  background  signal 
for  an  ideal,  imple  pulse  radar.  There  are  four  generic  types  of  back¬ 
ground  which  can  be  simulated.  The  first  is  white,  Gaussian  receiver 
noise,  which  is  produced  by  a  call  to  subroutine  NOISE.  The  second  and 
third  are  generalizations  of  NOISE,  wherein  the  amplitude  statistics 
of  the  receiver  noise  are  changed  from  Rayleigh  to  either  log-normal 
(LGNORM)  or  Weibull  (WEIBUL) .  In  all  cases ,  samples  of  backgrounds 
produced  by  these  three  generators  are  independent  from  sample  to  sample. 
The  fourth  type  of  background  generator,  CORBKG, is  an  implementation 
of  the  clutter  model  of  section  2.2.  As  such,  it  generates  random 
samples  of  average  clutter  cross-section  with  either  log-normal  or 
Weibul  statistics,  and  then  applies  uncorrelated  Rayleigh  amplitude 
and  uniform  phase  modulations  to  each  range  bin  from  pulse- to-pulse. 


31 


DETDFX 


PHCODE 

TARGET 

PDDFX 

SUMARY 

I 

1  ZZ1 _ 

INIT 

DBCONV 

OPNCLS 

PLOT 


Figure  18  -  Top  Level  Simulation  Organization 
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The  waveform  and  simple  pulse  radar  backgrounds  are  combined 
Co  form  a  bipolar  video  signal  for  a  pulse  compression  radar  with  a  call 
Co  IFGEN.  This  subroutine  convolves  Che  transmit  waveform  produced  by 
PNCODE  with  the  backgrounds  from  the  generators.  The  convolution  -is 
in  the  range  dimension  only.  When  the  background  is  Gaussian  noise, 
the  call  to  IFGEN  is  by-passed. 

Detection  statistics  are  accumulated  and  summarized  by  the  three 
modules  TARGET,  PDDFX,  and  SUMARY.  TARGET  is  a  target  generator  module, 
which  returns  samples  of  target  amplitude  which  are  either  non-fluctuating 
(Swerling  case  0),  or  one  of  the  four  Swerling  cases  1-4.  PDDFX  is 
the  module  which  implements  the  hard-limited  pulse  compression  oper¬ 
ations.  It  first  adds  the  target  signal,  assembled  from  the  waveform, 
the  target  amplitude  generator,  and  an  input  value  of  signal-to-background 
ratio,  to  the  linear  signals  generated  by  IFGEN.  The  sum  is  then  hard- 
limited  and  pulse  compressed.  In-phase  and  quadrature  channel  outputs 
are  squared  and  summed,  and  combined  with  the  corresponding  outputs  from 
other  diversity  channels.  The  resultant  is  then  subjected  to  a  threshold 
test  at  each  of  several  threshold  levels.  This  process  of  target  and 
background  addition,  pulse  compression,  diversity  channel  combination, 
and  threshold  comparison  is  repeated  for  all  target-to-background  ratios. 

A  matrix  of  threshold  crossing  counts  is  accumulated  by  repeating  for 
many  Independent  trials.  Module  SUMARY  summarizes  these  counts  in  a 
tabular  format,  expressed  as  probability  of  threshold  crossing,  as  a 
function  of  threshold  level  and  signal-to-background  ratio.  The  PLOTPD 
module  writes  the  threshold  crossing  probabilities  to  a  file  for  off-line 
plotting.  Additionally,  it  obtains  a  parametric  fit  to  the  data,  writes 
the  parameter  values  to  the  output  file,  and  generates  and  writes  to 
the  plot  file  the  smooth  fit  curve. 

The  three  remaining  modules  are  utilities  for  array  initialization 
(INIT) ,  linear  to  dB  conversion  (DBCONV) ,  and  output  file  opening  and 
closing  (OPNCLS). 
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4.2  Simulation  Simplifications 

As  in  any  simulation,  a  number  of  simplifying  assumptions  have 
been  invoked  in  the  interest  of  efficiency.  Some  of  the  more  significant 
simplifications  are  mentioned  here. 

4.2.1  Waveform  Selection 

As  previously  mentioned,  the  simulation  is  restricted  to  wave¬ 
forms  consisting  of  0°  and  180°  phase  code  modulations  of  unity  amplitude 
sinusoidal  sub-pulses.  The  phase  codes  are  generated  from  maximal  length 
sequences  with  shift  register  generators.  In  general,  there  are  several 
alternate  shift  register  feedback  connections  which  will  generate  a  given 
length  code.  In  addition,  for  each  such  set  of  connections,  there  is 
an  arbitrary  starting  point  in  the  resulting  maximal  length  sequence, 
corresponding  to  the  initial  state  of  the  shift  register  generator. 

The  classic  approach  to  selecting  from  among  all  such  possible 
resulting  waveforms  of  a  given  length,  is  to  examine  the  linear  processing 
ambiguity  function  (point  target  matched  filter  response  as  a  function 
of  delay  and  Doppler  mismatch)  of  each  waveform,  and  select  one  which 
exhibits  the  "best"  sldelobe  behavior.  In  the  non-linear  processing 
schemes  simulated  here,  this  is  still  probably  a  reasonably  good  selection 
criterion.  However,  in  the  present  simulation,  an  arbitrary  selection 
of  sequence  has  been  made  by  using  the  shift  register  feedback  con¬ 
nections  suggested  in  Table  6  of  Chapter  20  of  the  Radar  Handbook  [3], 
and  an  Initial  shift  register  state  of  all  ones.  The  autocorrelation 
function  sidelobes  of  this  waveform  are  plotted  in  Figure  19.  The  peak 
sldelobe  is  about  19  dB  below  the  mainlobe. 

4.2.2  In-phase  Signal  Assumption 

The  target  signal  return  is  always  assumed  to  have  an  initial 
RF  phase  of  0°,  so  that  all  target  energy  occurs  in  the  in-phase  arm 
of  the  I  and  Q  channel  processor.  This  differs  from  the  actual  situation 
in  which  target  phase  will  be  uniformly  distributed  over  360°,  with 
target  energy  split  between  the  two  quadrature  channels.  For  linear 
processing,  this  assumption  of  zero  target  phase  has  no  effect  on  the 
performance  of  the  unknown  phase  (I  and  Q  channel)  receiver,  since  the 
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receiver  Is  equivalent  to  matching  to  an  arbitrary  phase  followed  by 
envelope  detection  [4]. 


The  zero  target  phase  assumption  is  believed  to  have  little 
effect  on  performance  of  the  hard- limiting  receiver  also,  although 
this'  assumption  has  not  been  justified. 

4.2.3  Beamshape,  Range  Bin,  and  Doppler  Bin  Scalloping  Losses  Neglected 

A  target  seldom  falls  on  the  main  response  axis  of  the  radar 
beam,  nor  is  it  always  aligned  with  the  center  of  a  range  or  Doppler  bin. 

As  the  beam  sweeps  past  the  target,  the  sequence  of  target  returns  is 
modulated  by  the  beamshape,  with  the  result  that  the  single  pulse  slgnal- 
to-noise  ratio  is  not  constant  within  the  pulse  sequence.  A  similar 
effect  occurs  due  to  mismatch  in  range  and  Doppler,  often  called 
scalloping. 


These  effects  have  been  neglected  in  the  simulation.  Target 
amplitude  is  constant  throughout  the  sequence  of  N  pulses.  Similarly, 
the  target  is  presumed  to  fall  in  the  center  of  a  range  bin,  for  all 
pulses.  Finally,  target  Doppler  is  assumed  to  be  zero. 
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APPENDIX  A 


A.O  DETAILED  SIMULATION  PROGRAM  DESCRIPTION 

This  appendix  contains  a  more  detailed  description  of  the  simulation 
program  than  is  contained  in  the  main  body  of  the  text.  The  description  is 
done  in  two  parts.  The  first  part  is  module-by-module  description  of  the 
entire  program,  in  which  a  narrative  discussion  is  given  of  the  computation 
performed  by  each  piece  of  the  program.  This  is  essentially  an  elaboration 
of  the  functional  summary  contained  in  Section  4.1.  Following  this,  an 
example  run  is  presented  showing  inputs  and  outputs. 

A.l  Module  Description 

Figure  A.l  presents  the  module  hierarchy  of  the  simulation  program, 
showing  the  relationship  of  all  subroutines.  It  is  analagous  to  the  organi¬ 
zation  chart  shown  in  Figure  18,  but  contains  greater  detail,  and  presents 
the  subroutines  in  the  order  of  call.  Each  module  in  this  chart  is  described 
in  the  paragraphs  which  follow.  Section  A. 2  contains  the  example  run. 
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Figure  A.l  Module  Hierarchy 
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*  System  Routine 


A.  1.1 


Driver  (DETDFX) 


Calling  Statement:  driver  module 

Parameters  Read  from  Input  File*; 

NTRIAL  -  number  of  Monte-Carlo  range-bins  for  simulation  run 
NPDI  *  number  of  pulses  (or  frequency  diversity  channels) 
for  incoherent  integration 
M  ■  pulse  compression  ratio 

IPLOT  •  flag  to  specify  generation  of  a  plot  file  for  off¬ 
line  plotting  of  results 
IPLOT  *  1  ■  generate  plot  file 
IPLOT  +  1  -  no  plot  file 
NTRLPS  *  input  number  of  range  bins  per  scan. 

ITYPE  -  background  type  code 

ITYPE  ■  1  •  Gaussian  receiver  noise. 

ITYPE  «  2  ■  Log-normal  uncorrelated  clutter 
ITYPE  *  3  ■  Weibull  uncorrelated  clutter 
ITYPE  -  4  *  Correlated  clutter 

NTHR  ■  number  of  threshold  levels  for  accumulation  of 
threshold  crossing  counts 
(THR(I) ,  1-1,  NTHR)"  10  log^g  (threshold  levels) 

NSBR  "  number  of  signal-to-background  ratios 
(SBR(I),  1-1,  NSBR)"  signal-to-background  ratios  in  dB. 


Description : 

The  driver  module  (DETDFX)  performs  two  functions.  It  serves 
as  the  executive  routine  which  defines  the  sequence  of  functions  needed 
for  simulation  of  the  signal  processor  described  in  Section  2.1.  It  also 
performs  a  limited  amount  of  "book-keeping"  calculations  necessary  for 
proper  executive  control. 

The  first  operation  is  a  call  to  subroutine  OPNCLS,  which 
opens  input  and  output  files  on  Fortran  logical  units  LUIN  and  LDOUT. 

All  program  Inputs  are  made  on  LUIN.  With  the  exception  of  plot  file 
outputs  produced  by  subroutine  PLOTPD,  all  program  outputs  are  directed 
to  LUOUT. 
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input  parameter  types  and  ranges  of  allowable  values  are  summarized  in 
Section  A. 2. 
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An  unformatted  read  from  LUIN  seta  values  for  parameters  NTRIAL 
NPDI,  M,  IPLOT,  and  NTRLPS.  These  parameters  are  defined  briefly  above. 

A  value  of  zero  for  NTRIAL  signals  the  end  of  a  simulation  run. 

For  non-zero  NTRIAL,  the  next  step  is  a  call  to  subroutine 
IN1T,  which  Initializes  all  arrays  in  labelled  conmon  areas  CNTRL,  ARRAYS 
and  FLAGS.  The  total  number  of  range  bins  (NTRIAL)  will  be  simulated  in 
blocks  of  NRBIN  range  bins  at  a  time.  This  will  be  repeated  for  NSCAN 
blocks,  until  the  specified  total,  NTRIAL,  is  reached.  As  a  first  cut, 
these  two  factors  are  set  from  the  input  variables  NTRIAL  and  NTRLPS 
according  to: 

NSCAN  •  NTRIAL/NTRLPS 
NRBIN  -  NTRLPS 

This  initial  value  of  NRBIN  is  now  checked  to  see  if  it  requires  NBIN, 
the  number  of  complex  sample  pairs  per  call  to  the  background  generator 
NOISE,  LGNORM,  WEIBUL,  and  CORBKG,  greater  than  2048.  If  not,  than  these 

initial,  values  are  acceptable .  Otherwise  NRBIN  is  descremented  and  NSCAN 
recalculated: 

NSCAN  -  NTRIAL/NRBIN 

This  is  repeated  until  NBIN  is  less-than-or-equal-to  2048.  When  the 
constraint  on  NBIN  is  satisfied,  NSCAN  and  NRBIN  are  checked  to  see  if 
they  produce  at  least  the  desired  number  of  trials.  If  not,  NSCAN  is 
incremented,  and  NTRIAL  re-calculated. 

The  remaining  inputs,  ITYPE,  NTHR,  THR,  NSBR,  and  SBR  are  now 
read  in  from  LUIN,  and  a  summary  of  Inputs  is  written  to  LUOUT. 

A  call  to  PNCODE  generates  an  M-sample  representation  of  the 
transmitted  waveform,  returning  it  in  array  WR.  All  samples  of  signal 
returned  by  PNCODE  are  either  +1  or  -1. 


The  next  seep  is  a  conversion  from  dB  units  to  linear  units  of 
THR  and  SBR.  The  threshold  conversion  is  'straightforward,  according  to 

THR  -  10<IHR/1I» 

The  conversion  of  SBR  is  such  that  the  converted  value  becomes  the  target 
amplitude  which  would  produce  the  desired  signal-to-background  ratio  at  the 
output  of  a  linear  matched  filter  pulse  compressor.  This  conversion  depends 
on  the  type  of  background.  For  receiver  noise,  the  conversion  is 

..(SBR/ 20) 

SBR  -  — - 

/F 

For  distributed  clutter  backgrounds,  the  conversion  is  simply 

SBR  -  10(SBR/20) 

Following  conversion  of  SBR  and  THR,  the  main  loop  of  the  simulation 
begins.  For  each  scan,  a  matrix  of  background  samples  is  generated  with  a 
call  to  the  appropriate  background  generator.  The  samples  are  returned  in 
VUR  and  VUI  by  the  receiver  noise  generator,  and  in  XNR  and  XNI  by  the 
distributed  clutter  generators.  The  sample  matrix  is  large  enough  such  that 
NRBIN  pulse-compressed  range  bin  outputs  can  be  simulated  on  each  of  NPDI 
pulses  (or  diversity  channels) .  A  call  to  IFGEN  generates  a  simulated  radar 
bipolar  video  signal  by  convolving  the  background  samples  with  the  waveform 
in  the  range  dimension.  The  video  samples  are  returned  in  the  arrays 
VUR  and  VUI.  If  the  background  type  is  receiver  noise,  this  convolution 
step  is  by-passed.  A  call  to  TARGET  generates  NRBIN  x  NPDI  target  samples, 
which  are  returned  in  array  TGT. 

Subroutine  PDDFX  combines  target  and  background  signals,  hard- 

limits  the  in-phase  and  quadrature  resultant  signals,  and  pulse  compresses. 

This  is  repeated  on  each  range  bln  in  a  scan,  on  each  of  NPDI  pulses  (or 

diversity  channels) .  The  pulse  compressed  outputs  are  square-law  detected 
2  2 

(I  +  Q  )  and  a  test  statistic  generated  by  summing  over  pulses.  The  test 
statistic  is  then  compared  to  the  vector  of  threshold  levels  THR,  and 
crossing  counts  accumulated  in  NTC  for  each  threshold  and  signal-to-background 
ratio  combination. 


Following  the  major  loop,  THR  and  SBR  are  re-converted  to  dB,  and 
the  threshold  crossing  counts  are  summarized  as  probability  of  threshold 
crossing  (P^)  by  a  call  to  SUMARY.  An  optional  call  to  PLOTPD  produces  a 
parametric  fit  to  each  P^  vs  SBR  curve,  and  writes  the  simulated  data  points 
and  curve  fit  to  a  file  for  off-line  plotting. 

The  program  then  loops  back  to  the  beginning  for  another  case. 

When  all  cases  are  finished  (NTR1AL  *  0),  the  program  ends  with  a  final  call 
to  OPNCLS  to  close  LUIN  and  LUOUT. 
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C  PROGRAM  TO  SIMULATE  DETECTION  OF  TARGETS  IN 
C  NOISE  OR  CLUTTER  FOR  DICKE-FIX  CFAR. 

C 

INTEGER *4  NTRIAL, NSC AN 
* INSERT  JIMB  JOD I >SUR ADS>CNTROL 
-SINSERT  JIMB  JODI>SURADS>ARRAYS 
<5 INSERT  JIMB  JODI>SURADS>IOLUS 

C 

DATA  INEWPG/:  214/ 

CALL  OPNCLS(LUIN<  LUOUT*  1 ) 

C 

1  WRITE (LUOUT.  100) 

100  FORMAT < ////IX.  'DETECTION  PERFORMANCE  OF  DICKE-FIX  CFAR') 

C 

READ (LU IN. *)  NTRIAL. NPDI.  M» IPLOT. NTRLPS 
IF (NTRIAL.  EQ.  0)  GO  TO  40 
C 

CALL  INIT 

C 

MM1T2-(M-1>*2 

NSC AN-NTR I AL /NTRLPS 

NRB I N-NTRLPS 

5  CONTINUE 
NRM2-NRB I N+MM 1 T2 
NBIN-NRM2*NPDI 

IF < NB IN.  LE.  2048)  GO  TO  6 

NRBIN-NRBIN-1 

NSC  AN— NTR  IAL/NRBIN 

GO  TO  5 

6  CONTINUE 

IF( ( NSC AN*NR BIN).  LT.  NTRIAL)  NSCAN-NSCAN+1 
NTR I AL»NSCAN*NRB IN 

C 

READ(LUIN. *)  I TYPE.  NTHR,  (THR(I).  1-1,  NTHR). NSBR.  (SBR(I).  I-l. 
WR I TE ( LUOUT .  200)  NSCAN, NRB IN.  IPLOT,  I TYPE,  NPDI,  M 
200  FORMAT  ( ///6X,  'NSCANS  -  ',  15,  5X,  'NRB  I  NS  -  ',  15.  5X,  'IPLOT  -  ' 
*///6X,  'BACKGROUND  TYPE  »  ',  13,  5X,  'NPDI  »  ',  13,  3X,  'M  -  '.  13) 
WR I TE ( LUOUT ,  300)(THR(I),  I-l,  NTHR) 

300  FORMAT ( / /6X, 'THRESHOLD  TEST  LEVELS  '/ 

*5(/6X,  4011.  3/)  ) 

WR I TE ( LUOUT ,  400)  (SBR(I),  I-l,  NSBR) 

400  FORMAT ( //6X, 'SIGNAL-TO-BACKGROUND  RATIOS  (DB)  :  '/ 

*5(/6X, 4G11. 3/) ) 

C 

CALL  PNCODE(WR.M) 

C 

C  CONVERT  S/B  RATIOS  TO  LINEAR  UNITS 

C 

CALL  DBCONV(SBR, NSBR,  .  05,  1 ) 

CALL  DBC0NV1THR, NTHR,  .1,1) 

IF ( .  NOT.  (ITYPE.  EQ.  1 )  )  GO  TO  20 
SAMP-1.  /SORT (FLOAT (M) ) 

DO  10  I-l,  NSBR 

SBR ( I ) =SAMP*SBR ( I ) 

10  CONTINUE 

20  CONTINUE 

C 

DO  30  I-l, NSCAN 

IF< ITYPE.  EQ.  1 )  CALL  NOISE(VUR, VUI, NBIN, -1 ) 

IF( ITYPE.  EQ.  2)  CALL  LGNORMl XNR, XNI, NBIN, -1 ) 

IF( ITYPE.  EQ  3)  CALL  WEIBUL( XNR, XNI,  NBIN, -1 ) 

IF(  ITYPE.  EQ.  4)  CALL  CORBKG(  XNR,  XNI,  CR,  NRM2,  NPDI, -1,  1 ) 

IF ( ITYPE.  NE.  1 ) 

K  CALL  IFGEN  (VUR, VUI, XNR, XNI,  NRM2,  NPDI,  WR,  M) 

CALL  TARGET (TGT, NRB IN, NPDI ) 

CALL  PDDFX(M) 

30  CONTINUE 

C 

C  CONVERT  S/B  RATIOS  TO  DB 

C 

IF < .  NOT.  ( ITYPE.  EQ.  1 )  )  GO  TO  36 
DO  32  I-l, NSBR 
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UL 


NSBR) 
,  15, 


SBR ( I >  »SBR ( I ) /SAMP 
32  CONTINUE 

36  CONTINUE 

CALL  DBCONV<  SBR, NSBR,  20.  i  -1 ) 

CALL  DBCONVtTHR, NTHR,  10.  ,-l) 

C 

WR ITE ( LUOUT, 460 )  INEWPG 
460  FORMAT <A2> 

CALL  SUMARY <NTC, SBR,  THR, NSBR,  NTHR, NTRIAL) 

IF< I PLOT.  EQ.  1 )  CALL  PLOTPD<NTC,  SBR,  NSBR,  NTHR,  NTRIAL) 
WR I TE  <  LUOUT , 300 )  I NEWPO 
300  FORMAT (A2/70 <'*') ) 

eo  to  l 

c 

40  CALL  OPNCLS<LUIN,  LUOUT,  0) 

STOP 


*  * 


INTEGER *4  NTC* NSC AN 

COMMON/CNTROL/NTC  (  400  > .  SBR  <  20  > .  THR-<  20  > , 
*  NSBR.NTHR.NSCAN,  NRBIN,  NPDI.CNR 


COMMON/ ARRAYS/CR < 2048 ) ,  C I < 2048  >  * XNR <  2048 ) »  XNX  < 2048 ) » 

VUR  < 2048 )  i  VUI (2048) .  VCR (2048 ) » VC I <2048) < 
TOT (2048).  WR(2S6>< MI (296) 


COMMON/ IOLUS/LUIN,  LUOUT 
C 


COMMON/FLAGS/ IFLAON<  IFLAGL* IFL.AOW,  IFLAOT,  IFLAO(IOO) 
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A. 1.2 


File  Open  and  Close  (OPNCLS) 

Calling  Statements:  Call  OPNCLS  (LOIN,  LUOUT,  ICODE) 


Argument  List  Definition: 

Input  Variables:  ICODE  ■  open/close  code 

(1  ■+•  open,  0  close) 

Returned  Variables:  LUIN  -  input  logical  unit 

LUOUT  *  output  logical  unit 

Other  Modules  Called:  SEARCH,  TIMDAT 

Description: 

OPNCLS  opens/closes  input/output  files,  using  system  subroutine 
SEARCH,  if  job  is  a  batch  job  (NUSER  10).  Input  file  is  'CFARIN',  output 
file  is  'CFROUT' .  LUIN  and  LUOUT  are  set  to  5  and  6,  respectively. 

If  job  is  conversational,  LUIN  and  LUOUT  are  set  to  1  and  2,  and 
all  I/O  is  to  user  terminal. 

In  addition,  OPNCLS  obtains  time  and  date  of  run  from  system 
subroutine  TIMDAT,  and  documents  on  LUOUT. 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

10 

c 

c 

c 


c 

20 

c 

c 

c 


900 


c 

30 

C 

c 

c 

c 

c 

c 


c 

40 

100 

200 

300 

400 


SUBROUTINE  OPNCLS( LUIN,  LUOUT,  I CODE) 

SUBROUTINE  TO  SE1  FORTRAN  LOGICAL  UNITS 
LUIN  AND  LUOUT.  AND  OPEN  AND  CLOSE  DISK 
INPUT  AND  OUTPUT  FILES  FOR  PHANTOM  JOBS. 

INTEGER  OUTFIL,  OU I DSK 

DIMENSION  INFO( IS),  INFILE(3) .  0UTFIL<3> 

DATA  INFILE.  OUTFIL/ 'CF',  'AR',  'IN'.  'CF'.  'RO'.  'UT'/ 
DATA  INDISK.  OUTDSK/1.  2/ 

DATA  NPHAN/10/ 

IF(.  NOT.  ( I  CODE.  EQ.  1 )  >  GO  TO  30 

OPEN  FILES 

CALL  TIMDATdNPO,  13) 

NUSER-INF0(12> 

IMIN— MOD(  INFQ(4) »  60) 

I HRS- INFO (4) /60 
ISEC— INFQ(S) 

IF < .  NOT.  (NUSEN.  GE.  NPHAN) )  GO  TO  10 

PHANTOM  JOB: 

LUIN— INDISK+4 
LUOUT -0UTD8K+4 

CALL  SEARCHd,  INFILE.  INDISK.  *100) 

CALL  SEARCH <2.  OUTFIL. OUTDSK.  4200) 

GO  TO  20 

CONTINUE 

INTERACTIVE  JOB: 

LUIN— 1 
LUOUT— 1 

CONTINUE 

WRITE  TIME  AND  DATE  OF  RUN 


WRITE<LU0U1. 900)  (INFO(I),  1-1,3), 

*  (INFO(I). 1-13,  IS),  INF0(12) 

FORMAT ( ////IX,  S0(  '* ' )  /SX,  A2,  2(  '/'< 

*  2X.3A2,  '  (USER  # ',  12,  '  > '/IX,  S0( '* 
GO  TO  40 

CONTINUE 

CLOSE  FILES 

IF ( .  NOT.  (NUSER.  GE.  NPHAN) )  GO  TO  40 
PHANTOM  JOB: 

CALL  SEARCH ( 4,  INFILE,  INDISK,  4300) 
CALL  SEARCH( 4, OUTFIL,  OUTDSK,  *400) 

CONTINUE 
RETURN 
STOP  1 
STOP  2 
STOP  3 
STOP  4 
END 


I  HRS.  IMIN,  ISEC 

,  A2) , 2X,  2(12,  ' 
')/) 
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A. 1.2.1  System  Routine  TXMDAT 

********** 

*  TXMDAT  * 

********** 


The  calling  sequence  is: 

CALL  TXMDAT  (array,  nun) 

TXMDAT  returns  the  date,  time,  CPU  tine,  and  paging  tine  used  since 
LOGIN,  the  users  unique  nunber  on  the  system,  and  his  login  UFD  name  in 
an  array  as  follows: 

array  (1)  Two  ASCII  characters  representing  month.  Example:  11 

(2)  TWo  ASCII  characters  representing  day.  Example:  30 

(3)  TWo  ASCII  characters  representing  year 
Example:  75 

(4)  Integer  time  in  minutes  since  midnight. 

(5)  Integer  time  in  seconds. 

(6)  Integer  tine  in  ticks. 

(7)  Integer  CPU  time  used  in  seconds. 

(8)  Integer  CPU  tine  used  in  ticks. 

(9)  Integer  disk  I/O  time  used  in  seconds. 

(10)  Integer  disk  I/O  tine  used  in  ticks. 

(11)  Integer  nunber  oL  ticks  per  second. 

(12)  User  nunber. 

(13)  Six-character  login  name,  left-justified. 

(14) 

(15)  Example:  MSMITB 


num  words  of  array  are  set 


A.  1.2. 2  System  Routine  SEARCH 


SEARCH 


Definition  of  SEARCH 

SEARCH  is  used  to  connect  a  file  to  a  file  unit  (cpen  a  file)  or 
disconnect  a  file  from  a  file  unit  (close  a  file) .  After  a  file  is 
oormected  to  a  unit,  FEWFIL  and  other  routines  may  be  called,  either  to 
position  the  current-position  pointer  of  a  file  unit  (file  pointer)  or 
to  transfer  information  to  or  from  the  file  (using  the  file  unit  to 
reference  the  file) . 

Opening  a  File 

Ctt  opening  a  file,  SEARCH  specifies  1)  allowable  operations  that  may  be 
performed  by  FBHFXL  and  other  routines  (these  operations  are  read  only, 
write  only,  or  both  read  and  write) ;  2)  where  to  look  for  a  file  or 
where  to  add  the  file,  if  the  file  does  not  already  exist;  and  3) 
whether  the  file  is  to  be  opened  for  writing  only  or  both  reading  and 
writing.  SEARCH  either  specifies  a  filename  in  the  currently  attached 
user  file  directory  or  a  file  unit  number  on  which  a  segment  directory 
is  open.  In  the  segment  directory  reference,  the  file  to  be  opened  or 
closed  has  its  beginning  disk  address  given  by  the  word  at  the  current 
position  pointer  of  the  file  unit. 

SEARCH  Actions 

Cn  creating  a  new  file,  the  user  specifies  to  SEARCH  the  file  type  of 
the  new  file. 

The  subroutine  SEARCH  may  be  used  to  perform  actions  other  than  opening 
and  closing  a  file.  SEARCH  may  delete  a  file,  rewind  a  file  unit,  or 
truncate  a  file. 

Upon  opening  a  file,  SEARCH  sets  the  file  pointer  to  the  beginning  of 
the  file.  Subroutines  P5WFIL  and  others  cause  information  to  be 
transfered  to  or  from  the  file  unit,  starting  at  the  file  pointer. 
After  the  transfer,  the  pointer  is  moved  past  the  data  transferred.  A 
call  to  SEARCH  to  rewind  a  file  causes  the  file  pointer  to  be  set  to 
the  beginning  of  the  file.  Subsequent  calls  to  FEWFIL  and  other 
routines  cause  information  transfer  to  occur  as  if  the  file  had  just 
been  opened.  A  call  to  SEARCH  to  truncate  a  file  causes  all 
information  beyond  the  file  pointer  to  be  removed  from  the  file.  This 
call  is  useful  if  one  is  overwriting  a  file  with  less  information  than 
was  originally  contained  in  the  file. 
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Subroutine  Call 

SEARCH  is  used  as  in  the  following  call: 

Format: 

CALL  SEARCH  (SSI,  NAME,  FUNIT,  AITKIN) 

KEY  is  composed  of  three  subkeys  that  are  combined  additively:  ACTION, 
REFERENCE,  and  NEWFILE.  Not  all  subkeys  are  required  on  every  call, 
and  subkeys  with  values  of  zero  can  be  omitted.  The  SEARCH  call  may 
therefore  be  represented  as: 


CALL  SEARCH  (ACTIOI+REFERENCE4NEWFILE,  NAME,  FUNIT,  AITKIN) 


All  require  an  ACTION  subkey.  The  ACTION  subkeys  are  shown  in 

the  following  table: 


ACTION  Octal  Value  Meaning 


OPNRED  1 

orawRT  2 

OJNBTH  3 

CLOSE  4 

ngrjTPB  5 

EXIST  6 

REWIND  7 

TSNCAT  IS 

QBACC  1000 


Open  NAME  for  reading  on  FUNIT 
Open  NAME  for  writing  on  FUNIT 
Open  NAME  for  both  reading  and  writing 
on  FUNIT 

Close  file  by  NAME  or  by  FUNIT 

Delete  file  NAME 

Check  to  see  if  file  exists. 

Rewind  file  on  FUNIT 
Truncate  file  on  FUNIT 
Change  access  of  file  to  FUNIT 


The  REFERENCE  subkeys  are  shown  in  the  following  table: 


REFERENCE  Octal  Value  Meaning 

UFDREF  0  Searches  for  file  NAME  in  the  current 

user  file  directory  (UFD)  (as  defined  by  a 
previous  ATTACH)  and  perform  the  action 
in  the  ACTION  subkey  on  the  specified  file. 

SEGREF  100  Performs  the  action  specified  in  the  ACTION 

subkey  on  the  file  with  the  location 
indicated  by  the  file  pointer  designated 
within  the  array  NAME  (1) . 

This  file  unit  must  be  an  open 
segment  directory. 

Only  those  calls  to  SEARCH  that  reference  a  file  in  a  UFD  or  segment 
director/  need  the  reference  key.  Calls  that  reference  file  units 
do  not  need  this  key. 
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The  following  table  lists  the  NEWPI L  subkeys: 


NEWPIL  Octal  Value  Meaning 


ttTFXLE  0 
NDFILE  2000 
NTSEG  4000 
ilDSEG  6000 
NEWUFD  10000 


Mew  threaded  (SAM)  file 
Maw  directed  (DAM)  file 
Mew  threaded  (SAM)  segment  directory 
Mew  directed  (DAM)  segment  directory 
Mew  Qser  File  Directory  (SAM) 


Only  those  calls  to  SEARCH  that  generate  a  new  file  require  a  NEWFIL 
subkey.  Ch  other  calls,  this  subkey  is  ignored. 


The  none  of  the  remaining  parameters  in  a  call  to  SEARCH  are  as 
follows: 

4 

NAME  If  the  reference  subkey  is  (JFDREF,  NAME  is  either  a  . 

six-character  Hollerith  expression  or  the  name  of  a  three" -word 
array  that  specifies  a  filename  (existing  or  not) . 

If  the  reference  subkey  is  OFDREF  and  NAME(l)  is  -1,  the 
current  OFD  is  opened.  NAME  ■  -1  must  be  used  only  in 
configuration  with  ACTION  subkeys  1,  2,  or  3.  Oner  status 
of  the  current  OFD  is  required. 

If  the  reference  subkey  is  SEG8EF,  NAME  is  a  file  unit  (1-16; 
'1-15  under  FRIMCS  II)  on  which  a  segment  directory  is  open. 

Ch  calls  in  which  the  ACTION  key  requires  only  a  file  unit  to 
specify  the  file  to  be  acted  on,  NAME  is  ignored  and,  usually, 
specified  as  0. 

FUN IT  Cn  calls  that  require  a  file  unit  number,  FUNTT  is  a  number 

1  to  16  (1-15  under  FRIMCS  II).  Ch  calls  that  require  no  unit 
number,  FUN IT  is  ignored  and  usually  specified  as  1. 

AITKIN  AITKIN  is  an  integer  variable  assigned  the  value  of  a  label 
return  in  the  user's  FORTRAN  program  to  be  used  as  an  alter¬ 
nate  in  case  of  uncorrectable  errors  (e.g.,  attempting  to 
open  a  file  that  is  already  open) .  If  this  argument  is  0  or 
omitted,  an  error  message  is  printed;  control  returns  to 
FRIMCS  if  any  error  should  occur  while  using  SEARCH. 


A. 1.3  Array  Initialization  (INIT) 

Calling  Statement:  Call  INIT 
Description: 

Initializes  signal  arrays  and  threshold  crossing  accumulators . 
Also  sets  flags  to  initialize  inverse  probability  tables  in  NOISE,  WEIBUL, 
LGNORM,  TARGET,  and  CORBKG. 
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SUBROUTINE  INIT 


A.  1.4 


Waveform  Generation  (PNCODE) 

Calling  Statement:  CALL  PNCODE  (WR,  M) 


Argument  List  Definitions; 

Input  variable:  M  ■  length  of  pseudo-noise  code 

to  be  generated 

Returned  variables:  WR  ■  vector  of  M  waveform  samples, 

all  of  which  are  either  +1  or  -1. 

Other  Modules  Called:  None 

Description 

PNCODE  generates  a  pseudo-noise  pulse  compression  waveform  derived 
from  a  maximal  length  shift  register  generator.  Such  generators  are  described 
in  Chapter  20  of  the  Radar  Handbook  [3].  The  length  M  of  the  waveform  (or, 
equivalently,  the  pulse-compression  ratio)  is  an  input  to  the  subroutine. 

This  length  is  then  converted  into  the  number  of  stages  in  the  shift  register 
using 

NSTAGE  -  log ^  (M  +  1) 

If  M  +  1  is  not  Initially  a  power  of  two,  it  is  rounded  down  to  the  next 
lower  such  power.  The  shift  register  is  then  initialized  to  all  one's. 

A  sequence  of  one's  and  zero's  is  generated  by  using  feedback  connections 
suggested  in  Table  6  of  Chapter  20  of  [3].  The  one’s  and  zero's  are  converted 
to  +1  or  -1,  respectively,  which  are  returned  as  waveform  samples.  Prior  to 
returning  to  the  calling  program,  the  shift  register  generator  and  resulting 
waveform  are  summarized  on  LUOUT. 


SUBROUTINE  PNCODE  (W,  M) 


GENERATES  A  PSEUDO-RANDOM  CODE  OF 
MAXIMAL  LENGTH  M  IN  W.  SEE  RADAR  HANDBOOK 
PAOE  20-20. 


•INSERT  JIMB  JOD I >SUR ADS> I OLUS 

DIMENSION  W(l>,  ISHREG(IO),  IC0NN<4,  10) ,  NCONN(  10) 

DATA  NCONN/O,  6*2,  4.  2.  2/ 

DATA  I CONN/O.  0.  0.  0.  2.  1.  0.  0.  3.  2.  0.  0.  4.  3.  0.  0.  9.  3.  0.  0. 

*  6.  3i  0.  0.  7.  6.  0.  0.  8>  6.  9.  4.  9.  9.  0.  0.  10.  7.  0.  0/ 

NSTAQE-ALOQ ( FLOAT < M+l ) )  /ALOG ( 2.  0)*.  01 
M— (2**NSTAGE)— 1 
NCON-NCONN <  NSTAGE ) 

DO  10  1-1 , NSTAGE 
ISHREG(I)-1 
10  CONTINUE 

ISHREGC NSTAGE )-l 
C 

DO  40  J-l.M 
C 

C  FORM  MOD  2  SUM  OF  FEEDBACK  STAGES 

C 

I SUM-0 

DO  20  K-l.NCON 

ISUM— ISUM+ISHREG( IC0M4(K.  NSTAGE) ) 

20  CONTINUE 

I SUM— MOD < I SUM.  2) 

C 

C  EXTRACT  OUTPUT 

C 

IOUT— ISHREG <  NSTAGE ) 

IF  (IOUT.  EG.  1)  W<J)-1.0 
IF  (IOUT.  EG.  0)  W(J>— 1.0 
C 

C  SHIFT  REOISTER  CONTENTS 

C 

I ST AGE— NST AGE 
DO  30  1-2.  NSTAGE 

ISHREG( ISTA6E)-ISHREG( 1ST AGE-1 ) 

I STAGE- I STAGE— 1 
30  CONTINUE 

ISHREG(1)-ISUM 
40  CONTINUE 

C 

C  WRITE  OUT  WAVEFORM 

C 

WRITE  (LUOUT. 100)  M.  NSTAGE. ( ICONN(K.  NSTAGE),  K— 1< NCON) 

100  FORMAT  (//IX, 'MAXIMAL  LENGTH  PN  SEQUENCE'/ 

*/6X,  'CODE  LENGTH  -  ',19/ 

*/6X,  'SHIFT  REGISTER  LENGTH  -  ',  19/ 

*/6X. 'FEEDBACK  CONNECTIONS  -  '.419) 


419) 


WRITE  (LUOUT.  200)  (W(J).U-l.M) 

FORMAT ( /6X.  'WAVEFORM:  ' /.  64( /6X.  16(  IX.  F3.  0) / ) ) 

RETURN 

END 


A.1.5 


DB  Conversion  (DBCONV) 

Calling  Statement;  CALL  DBCONV  (X,N,P, ICODE) 

Argument  List  Definitions: 

X  -  array  of  values  to  be  converted 
N  ■  length  of  X  array 

P  ■  power  law  to  use  In  conversion  (see  below) 

ICODE  «  direction  of  conversion 

ICODE  ■  1  ■  DB  to  linear  conversion 
ICODE  +  1  •  linear  to  DB  conversion 

Other  Subroutines  Called:  None 

Description; 

Converts  to  or  from  dB,  according  to  the  following  expressions: 
If  ICODE  -  1: 

X*-  10X‘P 

If  ICODE  +  Is 

X'-  P  •  log10(X) 
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SUBROUTINE  DBCONV< X.  N. P, I CODE) 
CONVERSION  TO  AND  FROM  OB 
DIMENSION  X(l) 

IF < .  NOT.  (ICOOE.  EQ.  1>>  GO  TO  20 

DO  10  I-t.N 

xm-io.  **(x<i>*p) 

CONTINUE 
GO  TO  40 

CONTINUE 
DO  30  I-l.N 

X ( I )«P*AL0C10( X ( I ) ) 

CONTINUE 

CONTINUE 

RETURN 

END 
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Gaussian  Noise  Generator  (NOISE  and  SETUPN) 

Calling  Statement:  CALL  NOISE  (XR,  XI,  N,  IREP) 

Parameters  Read  from  Input  File: 

ISEEDN  »  starting  seed  for  pseudo-random  number  generator 
Argument  List  Definition: 

Input  Variables:  N  »  number  of  complex  samples  generated 

by  each  call 

IREP  ■  representation  code 

Returned  Variables:  XR  and  XI  are  defined  according  to 

IREP,  as  discussed  below. 

Other  Modules  Called:  SETUPN,  RANAMP,  RANPHS,  WRTPIV 
Description: 

For  IREP  *  -1  NOISE  generates  N  complex  samples  of  Gaussian  random 
noise  of  unit  mean  power,  and  uniformly  distributed  phase.  That  is, 

XR  -  /T  sin  $ 

XI  -  /F  cos  $ 

where  Y  is  exponentially  distributed  with  probability  density 

p(Y)dY  -  exp{-Y>  ,  Y  >  0 

and 

<J>  -  tan  ^(XR,  XI)  »  uniform  on  (-ir,ir). 

E[XR2  +  XI2]  -  E[Yj  -  1.0 

It  follows  that  XR  and  XI  are  uncorrelated  samples  from  a  Gaussian  distribution, 
with  density 

l  -t2 

p(X)dx  - -  exp{-=x) 

/2n  a  2  a 


where 


and 


X  is  XR  or  XI 

02  -  1/2. 
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In  IREP  *  -2,  NOISE  generates  sine  and  cosine  of  uniformly  distributed 
phase.  That  is, 

XR  -  sin  <j> 

XI  *  cos  $ 

For  IREP  “  -3,  XR  and  XI  are  hard-limited: 

XR  *  sign  (sin  $) 

XI  *  sign  (cos  <J>) 

For  non-negative  values  of  IREP,  NOISE  returns  detected  values  for 
square- law,  linear-lav,  and  logarithmic  detectors: 

IMP  ■  +2: 

irep-+1!  (Slf 

IMP  -  0:  I  Jn/r 

On  the  first  call  to  NOISE,  a  tabulation  of  the  inverse  cumulative 
distribution  for  /y",  Y,  or  ln/Y~  is  generated  and  stored  in  the  array  PIV. 

This  is  done  with  a  call  to  SETUPN.  The  table  is  a  three-level  table,  the 
first  level  evaluated  at  the  center  of  128  uniform  increments  of  the  0  to  1 
probability  axis.  The  second  level  divides  the  last  (nearest  unity)  such 
increment  into  128  sub-increments,  and  the  third  level  divides  the  last 
sub-increment  into  128  sub-sub-increments.  Thus  the  table  is  truncated  at 

probabilities  of  -tIf  *  .0039  and  1 - *3— ,  -  1  -  2.38  x  10  2 . 

128  (128)3 

Once  the  table  is  generated,  random  amplitude  samples  are  drawn 
from  the  tabulated  distribution  by  a  call  to  RANAMP.  Subroutine  RANPHS 
converts  to  Inphase  and  quadrature  components.  Subroutine  URTPIV  writes 
out  the  PIV  table  for  debugging  purposes. 


SUBROUTINE  NOISE (XR,  XI.  N,  IREP) 

C 

C  GAUSSIAN  NOISE  GENERATOR 

C 

C  IREP  RESULT 

C 

C  0:  RETURNS  NATURAL  LOGARITHM  OF  RAYLEIGH  AMPLITUDES  IN  XR. 

C 

C  1:  RETURNS  N  UNCORRELATED  RAYLEIGH  AMPLITUDES  IN  XR. 

C  OF  UNIT  MEAN  POWER. 

C 

C  2:  RETURNS  N  UNCORRELATED  EXPONENTIAL  VARIABLES  IN  XR. 

C  OF  UNIT  MEAN  POWER. 

C 

C  -1:  RETURNS  N  UNCORRELATED  COMPLEX  GAUSSIAN  SAMPLE  PAIRS. 

C  MEAN-0.  VARIANCE-O.  9.  RETURNS  N  SAMPLES  IN  XR. 

C  N  SAMPLES  IN  XI.  USES  3-LEVEL  TABLE  LOOK-UP  FOR 

C  AMPLITUDES. UNIT  CIRCLE  ALGORITHM  FOR  SIN  fc  COS./ 

C 

C  -2:  RETURNS  SIN  &  COS  (XR  &  XI >  OF  N  UNCORRELATED 
C  UNIFORMLY  DISTRIBUTED  PHASES. 

C 

C  -3:  RETURNS  SGN(SIN)  &  SGN(COS)  OF  N  UNCORRELATED  PHASES. 

C 

INTEGER *4  ISEEDN,  I SEED 
•INSERT  JIMB  AN I T A>SUR ADS> I OLUS 
•INSERT  JIMB  AN I T A>SUR ADS>FLAGS 

DIMENSION  XR(1), XI ( 1 ) , PIV(384) 

DATA  PIV.  LPIV. NPIV.  ISEED/384*1. 0. 3, 7.  1739483429/ 

C 

C  ON  FIRST  CALL  INITIALIZE  TABLE 

C 

IF  (.NOT.  (IFLAGN.  EG.  O)  )  GO  TO  29 
C 

WR I TE ( LUOUT .  100) 

100  FORMAT (////IX.  'GAUSSIAN  NOISE  GENERATOR') 

READ(LUIN. *)  ISEEDN 
IF ( ISEEDN.  EQ.  O)  I SEEDN— I SEED 
WRITE (LUOUT.  200)  ISEEDN 
200  FORMAT (//6X,  'PRN  GENERATOR  SEED  -  ',112) 

WRITE (LUOUT.  300)  IREP 

300  F0RMAT(//6X,  'REPRESENTATION  CODE-', 13) 

C 

CALL  SETUPN(PIV, NPIV,  LPIV,  IREP) 

IFLAGN— 1 
C 

29  CONTINUE 
C 

IF( IREP.  LT.  -1 )  GO  TO  28 
C 

C  GENERATE  RANDOM  AMPLITUDES  (OR  POWERS) 

C 

CALL  RANAMP(PIV,  NPIV.  LPIV,  ISEEDN,  XR,  N) 

C 

C  IF  IREP<0  GENERATE  AND  APPLY  RANDOM  PHASE 

C 

28  CONTINUE 

IF( IREP.  OE.  0)  GO  TO  30 

CALL  RANPHS( XR, XI,  N,  ISEEDN,  IREP) 

C 

30  RETURN 
END 


non 


SUBROUTINE  SETUPNtPIV,  NPIV,  LPIV,  IREP) 

SETUP  INVERSE  PROBABILITY  TABLE  FOR  GAUSSIAN  NOISE  GENERATOR 

DOUBLE  PRECISION  PROB 
•INSERT  JIMB  AN I T A>SUR ADS> I OLUS 
DIMENSION  PIV<1> 

C 

EXPO-.  5 

IF(  IREP.  EG.  2)  EXPO— 1.  0 

NN-2**NPIV 

JLP-0 

DO  20  LP-l.LPIV 
DO  10  J-l.NN 
JLP-JLP+1 

PROB- 1 .  DO+DBLE ( FLOAT ( J-NN >-.»>/<  DBLE ( FLOAT ( NN ) > **LP  > 
PIV( JLP  > -SNOL ( -DLOG ( 1.  DO-PROB >  > **EXPO 
IF ( IREP.  EG.  0)  PIV< JLP)»ALOG(PIV(JLP> ) 

10  CONTINUE 

20  CONTINUE 

C 

CALL  WRTPIV<PIV,  NPIV*  LPIV#  LUOUT) 

C 

30  RETURN 

END 


A.  1.7  PIV  Documentation  (WRTPIV) 


Calling  Statement:  Call  WRTPIV  (PIV,  NPIV,  LPIV,  LUOOT) 

Argument  List  Definition; 

Input  Variables:  PIV  ■  table  to  be  documented 

LPIV  «  number  of  levels  in  table 

NPIV  -  log2  (number  of  increments 
per  level) 

LUOUT  -  output  logical  unit 
documentation 


Description: 

A  call  to  WRTPIV  writes  out  the  contents  of  PIV  to  LUOUT,  using 
formatted  writes.  Useful  for  debugging. 


noon 


SUBROUTINE  WRTPIVtPIV,  NPIV.  LPIV,  LUOUT) 


C 


10 

C 

100 


20 

200 

30 

C 


WRITE  OUT  CONTENTS  OF  PIV. 

THE  INVERSE  PROBABILITY  TABLE. 

DOUBLE  PRECISION  P.DP.DP1.DP0 
DIMENSION  PIV<2),  P<9>,  DP<9),  X(9>.  ILABL(IO) 


NN-2**NPIV 
DP0-100.  DO 
DO  10  L-l.LPIV 

DP 1 -DPO/DBLE ( FLOAT ( NN  > ) 

P(L)»100.  D0-DP0-DP1*.  9D0 
DP(L)-DP1 
DP0-DP1 
CONTINUE 

WRITE (LUOUT, 100) < ( ILABLt I ) .  1-1,  10 > , L-l, LPIV) 
FORMAT ( ///6X,  'INVERSE  PROBABILITY  TABLE' 
*//6X, 9<10A2,  2X) // ) 

DO  30  N-l.NN 
INDX-N 

DO  20  L-l,  LPIV 
X(L>— PIV( INDX) 

P(L>«P(L)+DP(L) 

INDX-INDX+NN 

CONTINUE 

WR I TE ( LUOUT * 200) < <P(L),  X(L) ).  L-l,  LPIV) 
FORMAT  (AX,  9I0PF9.  9.  1P011.  3,  2X)  ) 

CONTINUE 

RETURN 

END 
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Random  Amplitude  Generator  (RANAMP) 

Calling  Statement:  Call  RANAMP  (PIV,  NPIV,  LPIV,  ISEED,  X,  K) 


Argument  List  Definition;  PIV  ■  inverse  cumulative  probability 

tabulation 

LPIV  ”  number  of  levels  in  table  (3  or  less) 

NPIV  ■  log. (number  of  increments  per 
level)  (NPIV  <  10) 

Input  Variables:  ISEED  *  pseudo  -  random  number  generator 

seed,  (double  length  Integer) 

K  -  number  of  samples  generated  per  call 

;  Output  Variables:  X  *  array  of  samples  from  tabulated 

distribution 

ISEED  *  modified  seed 


Other  Modules  Called:  NRAND,  SHFT,  INTS 


Description: 

RANAMP  implements  the  random  number  generation  via  multi-level 
table  look-up.  For  example,  for  a  three  level  table,  with  128  increments 
per  level,  random  amplitude  samples  are  drawn  from  the  distribution  tabulated 
in  PIV  by  generating  integers  uniformly  distributed  between  0  and  127,  and 
using  this  an  index  to  level  1  of  the  table.  When  a  value  of  0  is  drawn 
for  the  index,  a  second  random  Integer  is  drawn  to  index  level  2.  When 
the  second  ldex  comes  up  0,  a  third  integer  is  drawn  to  index  level  3. 

This  table  look-up  procedure  follows  one  by  Mitchell  [5] . 


For  each  sample  to  be  generated  a  32-bit  pseudo-random  integer  is 
generated  with  a  call  to  NRAND.  A  7-bit  positive  integer  is  obtained  from 
this  with  an  appropriate  shift-left,  shift-right  sequence.  If  this  integer 
is  non-zero,  it  is  used  as  the  index  to  a  level  1  table  look-up,  completing 
the  random  sample  generation.  If  it  is  zero,  7  more  non-over-lapping 
bits  are  shifted  out,  and  the  process  repeated  for  level  2  and  3  (if 
necessary)  of  the  table. 
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*v 


nnnnnooon 


C 


30 


35 

36 

40 

C 


SUBROUTINE  RANAHP(PIV,  NPIV.  LPIV, ISEEDi  X.K) 

SUBROUTINE  WHICH  GENERATES  RANDOM  VARIABLES 
USING  LPIV-LEVEL  TABLE  LOOK-UP.  PIV  IS 
TABLE  CONTAINING  INVERSE  CUMULATIVE 
PROBABILITY  DISTRIBUTION  FUNCTION.  TABLE 
LOOK-UP  ADDRESS  IS  GENERATED  WITH  CALL  TO 
PRN  SEQUENCE  GENERATOR  NR AND.  NPIV  BITS 
OF  PRN  ARE  USED  FOR  EACH  LOOK-UP. 

INTE0ER*4  I SEED,  NR AND. ILOC,  IL0C2 
DIMENSION  PIVUKXU) 

DATA  IBIT1/32/ 

NN-2**NPIV 
ISHRT-32-NPIV 
DO  40  I-l.K 

ILOC-NRAND(ISEFD) 

LEVEL- l 
IOFFST-O 
ISHLFT— IB IT1— 32 
CONTINUE 

IL0C2— SHFT ( ILOC,  ISHLFT,  ISHRT) 

IL0C3— INTS( IL0C2) 

IF  (ILOC  3.  NE.  0)  00  TO  36 

IF  (LEVEL.  EQ.  LPIV)  GO  TO  35 

LEVEL  -  LEVEL-*  1 

I OFFST - I OFFST +NN 

ISHLFT— ISHLFT-NPIV 

GO  TO  30 

CONTINUE 

IL0C3-NN 

CONTINUE 

I L0C3- I L0C3+ I OFFST 
X  (  I )  — P I V  ( IL0C3 ) 

CONTINUE 


RETURN 

END 
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A. 1.8.1  System  Function  SHFT 


SHFT  -  Logical  Shift  Function 

SHFT  performs  logical  shift  operations  on  integer  variables. 


Calling  Sequence: 

Format  1: 

IS  -  SHET  (I,  n»i) 

performs  a  shift  operation  on  the  variable.  If  IP1>0, 
tiie  shift  is  to  the  right;  if  IP1<0,  the  shift  is  too 
the  left;  if  IP  1-0,  no  shift  occurs.  This  operation 
is  equivalent  to  the  BS  function,  and  is  provided  far 
compatibility  with  other  FORTRAN  ccrtpilers. 

Format  2: 

IS  -  SHFT  (I,  IPl,  IP2) 

Format  2  performs  two  shift  operations,  first  by  IP1 
(setting  zeroes  in  vacated  bits) ,  then  by  IP2  (setting 
zeroes  in  vacated  bits) .  The  sign  of  IP2  and  IP2 
determine  the  direction  of  the  shift  while  their  mag¬ 
nitude  determines  the  number  of  bits  to  be  shifted. 


A.  1.8. 2  System  Function  INTS 

INTS  converts  its  argument  from  a  double  length  (32  bit)  to  a 
single  length  (16  bit)  Integer. 
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•  i. 


A.  1.9 


Pseudo-Random  Number  Generator  (NRAND) 

Calling  Statement:  ILOC  *  NRAND  (ISEED) 

Argument  List  Definition: 

Input  Variables:  ISEED  -  double  precision  integer 

pseudo-random  number  seed. 

Output  Variables:  NRAND  *  double  precision  Integer 

pseudo-random  number 

ISEED  -  NRAND 

Other  Modules  Called:  SHFT 
Description: 

NRAND  generates  a  pseudo-random  number  sequence  by  the  congruentlal 
method  [6].  That  is. 


X 


i+1 


■  [*  Vmod 


,32 


The  period  length  of  this  sequence  is  maximized  at 

29 

period  length  •  2  samples 


by  setting: 

x  odd 
o 

a  ■  8t  +  3,  t  -  1,  2,  3... 


The  value  for  t  which  has  been  used  here  is 

t-8192 


for  which 


a 


8  x  8192  +  3 


65539 

2 ^  +  2^  +  2 


0 


Hence,  the  multiplication  is  implemented  by  two  shifts  and  two  adds 


A. 1.10  Random  Phase  Generator  (RANPHS) 


Calling  Statement:  Call  RANPHS  UR,  XI,  N,  ISEED,  IREP) 

Argument  List  Definition; 

Input  Variables:  XR  -  input  vector  of  amplitudes 

(generated  by  RANAMP) 

N  “  Number  of  amplitude  samples 

ISEED  ■  Pseudo-random  No.  generator  seed 

IREP  -  representation  code 

Returned  Variables:  XR  and  XI  are  defined  according 

to  IREP,  as  discussed  below. 

ISEED  *  modified  seed 


Other  Modules  Called:  NRAND,  SHFT 
Description:  ’ 

RANPHS  generates  sine  and  cosine  of  a  uniformly  distributed  random 
phase  using  method  attributed  to  Von  Newman  by  Mithcell  [5].  The  procedure 
is  based  on  generating  (A,B)  coordinate  pairs  uniformly  distributed  within 
a  unit  circle,  and  computing  sine  and  cosine  terms  from  these  coordinates 
using  trlgonemtrlc  identities. 

The  coordinates  are  generated  from  a  32-bit  random  integer  obtained 

from  NRAND,  from  which  two  10-bit  2's  complement  integer  words  are  obtained. 

Each  such  word,  A  or  B,  is  uniformly  distributed  on  [-511,  512].  The  sum 

2 

of  the  squares  of  these  two  words  is  compared  to  (511)  -  261121,  to  test 
if  the  coordinate  pair  is  within  the  "unit"  circle,  where  "unit"  radius  is 
511.  If  so,  the  sine  and  cosine  terms  are  generated  using: 


cos  <p 


sin  <J> 


2AB 


2  2 
A*  +  o 


If  not,  another  coordinate  pair  is  generated  and  the  process  repeated 
until  a  point  in  the  unit  clrle  is  obtained. 


f 

■  ,r»  .  -■  ' 
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That  is: 


Tha  value  of 


detaralnas  the  returned  values  of  XR  and  XX 


IREP 

IREP  -  -1 

I REP  -  -2 

IREP  -  -3 


•  ( 


XR  > 
XI  ■ 


XR 

XI 


XR^sin* 
XR  C08$ 


sin$ 

cos$ 

sign(sin<f>) 

sign(cos4>) 


SUBROUTINE  RANPHS( XR, XI.  N, I SEED. IREP) 

GENERATES  SIN  AND  COS  OF  UNIFORM  PHASE. 

USING  UNIT  C1RCLF  ALGORITHM. 

IREP  RETURNS 
-1:  XI«XR*COS 

XR-XR*SIN 
-2:  XI-COS 

XR-SIN 

-3:  XI— *V-1 

XR-+/-1 

INTE0ER*4  I SEED.  INTOR. IRAND.  NR AND.  INTA.  INTB.  INTT.  INTT2 
DIMENSION  XR (  1  > .  XIU).  INTOR (2).  ISHL<2) 

EQUIVALENCE  t INTA.  INTGR ( 1 ) ) ,  ( INTB.  INTOR (2) > 

ISHRT-32-10 
ISHL( 1 ) -32-32 
ISHL<2>-22-32 

DO  40  I-l.N 

GENERATE  DOUBLE  PRECISION  PRN  INTEGER 
IRAND-NRAND < I SEED ) 

SHIFT  OUT  TWO  10-BIT  2'S  COMPLEMENT  WORDS 
(-912.  LE.  INTI.  LE.  911) 

DO  30  J-1,2 

ISHLF-ISHL(J) 

INTT-SHFT ( IRAND.  ISHLF) 

IF<  INTT.  0£.  O)  GO  TO  10 

INTT-NOT(INIT) 

.  INTT-SHFT ( INTT,  I SHRT ) 

INTT-NOT(INTT) 

GO  TO  20 

INTT-SHFT (INTT.  I SHRT) 

INTGR (J)«IN IT 
CONTINUE 

SQUARE.  SUM,  AND  TEST  UNIT  CIRCLE 
(911*911*261121) 

INTT«INTA*INTA 

INTT2*INTB«IN1B 

f SAND— INTT+INTT2 

IF( IRAND.  GT.  262121)  GO  TO  9 

IF  IREP— 3:  COMPUTE  SGN(SIN),  SGN(COS) 

IF ( .  NOT.  (IREP.  EQ. -3>)  00  TO  32 
XI ( I )— FLOAT ( ISIGN( 1 ,  INTA)) 

XR ( I )— FLOAT ( ISIGN( 1 » INTB)) 

GO  TO  40 

IF  IREP— 2:  COMPUTE  SIN  ,  COS 
CONTINUE 

IF(.NOT.  (IREP.  EQ. -2))  GO  TO  33 
R— 1.  0/FL0AT( IRAND) 

GO  TO  34 

IF  IREP— 1:  COMPUTE  X*8IN  AND  X*COS 
CONTINUE 

R-XR ( I ) /FLOAT ( IRAND ) 


34 

40 

C 


CONTINUE 

XI < I > -FLOAT  <2*INTA*INTB ) *R 
XR < X ) -FLOAT ( I N I T- I NTT2 > *R 
CONTINUE 


RETURN 

END 
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A.  1.11  Log-normal  Background  Ganerator  (LGNORM  and  SETUPL) 

Calling  Statement:  Call  LGNORM  (XR,  XI,  N,  IREP) 

Parameters  Read  from  Input  File; 

ISEEDL  ■  starting  seed  for  pseudo-random  number  generator 
SDB  -  log-std. -deviation  of  clutter  power  in  dB. 


Argument  List  Definition; 

Input  Variables:  N  -  number  of  complex  samples 

generated  by  each  call 

IREP  •  representation  code 

Returned  Variables:  XR  and  XI  are  defined  according  to 

IREP,  as  discussed  below. 


Other  Modules  Called:  SETXIPL,  RANAMP,  RANPHS 


Description : 

For  IREP  ■  -1,  LGNORM  generates  N  uncorrelated  complex  samples  of 
log-normal  video,  with  unit  mean  power  and  uniformly  distributed  phase. 

That  is: 

XR  -  /T  sin  <t> 

XI  -  /T  cos  $ 

where  Y  is  log-normally  distributed  with  probability  density 


p(Y)dY 


1 

S  y/2tF 

y 


-(InY  + 

exp  { - s- 

2S 

y 


,  Y  >_  0 


where 


S 

y 


E(Y) 

Y50 


log.*std.-dev.  of  Y 

/E [InY  -  InY]2 
SDB/4 .3429 

E(XR2  +  XI2)  -  1.0 

Sy2 

median  of  Y  -  exp  ( - *-) 
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and  where 


♦  •  tan  ^  (XR,  XI)  la  uniform  on  (-it, it). 

The  method  of  generation  of  the  complex  samples  Is  Identical  to  that  for 
NOISE.  An  Inverse  probability  table  Is  constructed  for  ft,  Y,  or  Inft, 
according  to  IRE?,  with  a  call  to  SETUPL.  A  sample  Is  drawn  from  this 
table  using  subroutine  RANAMP.  For  negative  IREP,  this  amplitude  sample 
is  then  converted  to  In-phase  and  quadrature  components  by  subroutine  RANPHS, 
which  generates  sine  and  cosine  of  a  uniformly  distributed  random  phase. 
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SUBROUTINE  LONQRM(XR, XI.  N.  IREP> 
LOG-NORMAL  CLUTTER  GENERATOR 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


RETURNS  NATURAL  LOGARITHM  OF  AMPLITUDES 
IF  IREP  WERE  1  (SEE  BELOW). 


1: 


RETURNS  N  UNCORRELATED  LOG-NORMAL  AMPLITUDES 
IN  XR,  WITH  UNIT  MEAN  POWER,  AND  LOO.  -STD.  -DEV. 
SDB( POWER  UNITS). 


RETURNS  N  UNCORRELATED  LOG-NORMAL  POWERS  IN  XR. 


-1: 


RETURNS  N  UNCORRELATED  LOG-NORMAL  CLUTTER  SAMPLES 
WITH  RANDOM  PHASE,  UNIT  MEAN  POWER,  AND 
LOG-STD. -DEV.  SDB.  RETURNS  NSAMPLES  IN  XR, 

AND  N  QUADRATURE  SAMPLES  IN  XI.  USES  3-LEVEL 
TABLE  LOOK-UP  FOR  AMPLITUDES,  UNIT  CIRCLE  ALGORITHM 
FOR  SIN  AND  COS. 


-2: 

-3: 


RETURNS  N  SIN  &  COS  PAIRS  ONLV,  IN  XR  8c  XI 
RETURNS  N  SGN(SIN)  8c  SGN(COS)  PAIRS  ONLY. 


4 INSERT  JIMB  AN I T A>SUR ADS> I OLUS 
♦INSERT  JIMB  AN I T A>SUR ADS>FL AGS 
INTEGER*4  ISEEDL,  I SEED 
DIMENSION  XR(1), XI(l), PI V (384) 

DATA  PIV, LPIV, NPIV, ISEED/384*!.  0,  3,  7,  298809973/ 


C 

C 

C 


C 

100 


ON  FIRST  CALL  READ  PARAMETERS  AND 
IF  (NOT.  (  IFLAGL.  EQ.  0)  )  GO  TO  23 


INITIALIZE  PIV  TABLE 


200 


WR I TE ( LUOUT ,  100) 

FORMAT (////IX, 'LOG-NORMAL  CLUTTER  GENERATOR 
READ(LUIN.  *)  ISEfr DL 
IF ( ISEEDL.  EQ.  O )  I SEEDL" I SEED 
WR ITE ( LUOUT,  200 )  ISEEDL 


300 


400 

C 


FORMAT ( //AX,  'PRN  GENERATOR  SEED  - 
READ(LUIN,  *)  SDB 
WRITE (LUOUT.  300)  SDB 

FORMAT (//AX,  'LOO.  -STD.  -DEV.  ■  ',  F6.  1 

WRITE (LUOUT.  400)  IREP 

FORMAT (//AX.  'REPRESENTATION  CODE  ■  ' 


112) 


'  DB') 
13) 


C 

23 


CALL  SETUPL(PIV, NPIV.  LPIV,  IREP,  SDB) 
IFLAGL"! 


C 

C 

C 


CONTINUE 

IF(  IREP.  LT.  -1) 


C 

C 

C 

28 


GO  TO  28 

GENERATE  RANDOM  AMPLITUDES  (OR  POWERS) 
CALL  RAN AMP (PIV, NPIV,  LPIV,  ISEEDL,  XR, N) 

IF  IREP-CO  GENERATE  AND  APPLY  RANDOM  PHASE 


C 

30 


CONTINUE 

IF ( IREP.  GE.  0)  GO  TO  30 

CALL  RANPHS( XR, XI, N,  ISEEDL,  IREP) 


RETURN 

END 
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non 


SUBROUTINE  SETUPLCPIV,  NPIV.  LPIV,  IREP.SDB) 

SETUP  INVERSE  PROBABILITY  TABLE  FOR  LOO-NORMAL  GENERATOR 


DOUBLE  PRECISION  PROB 
* INSERT  JIMB  AN I T A>SUR ADS> I OLUS 
DIMENSION  PIV(l) 

C 

EXPO-1.  0 

IF <  IREP.  EQ.  2>  EXPO-2.  0 
SZ-SDB/8.  6898 
NN— 2**NPIV 
JLP-0 

DO  20  LP-l.LPIV 
DO  10  J— 1>  NN 
JLP-JLP+1 

PROB— 1.  DO+DBLE < FLOAT ( J-NN ) -0 .  9) / < DBLEC FLOAT (NN) )**LP) 
P I V ( JLP ) -EXP ( EXPQ*SZ* <  OAUS I V ( PROB ) -SZ ) ) 

IF< IREP.  EQ.  0)  P I V <  JLP > - ALOO < P I V ( JLP ) ) 

10  CONTINUE 

20  CONTINUE 
C 

CALL  WRTPIV<PIV, NPIV,  LPIV,  LUOUT) 

C 

RETURN 

END 
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Inverse  Cumulative  Gaussian  Function  (GAUSIV) 

Calling  Statement:  X  ■  GAUSIV (P) 

Argument  List  Definition; 

Input  Variables:  P  *  probability  level  (double 

precision) 

Output  Variables:  GAUSIV  »  amplitude  for  probability 

level  P 

Description: 

This  function  is  a  rational  approximation  for  the  inverse  cummulative 
probability  distribution  function  for  zero-mean,  unit  variance  Gaussian  random 
variables.  It  is  based  on  equation  26.2.23  of  [7].  The  argument  P  must  be 
double  precision,  but  the  return  value  is  single  precision. 


A 
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FUNCTION  OAUSIVtZ) 

INVERSE  PDF  FOR  ZERO-MEAN-  UNIT  VARIANCE  9AUSSIAN  R.  V. 


DOUBLE  PRECISION  Z, X.  Y, C, D.  TEMP 
DIMENSION  C  <3> »  D<3) 

DATA  C/2.  919917,  .  802893,  .  010328/ 
DATA  D/1.  432788,  .  189269,  .  001308/ 
X«Z 


IF < Z.  OT.  0.  9D0)  X-l.DO-Z 
Y-DSQRT ( DLOG ( 1 .  DO/ ( X*X  >  > ) 
TEMP-C(1)+Y*(C<2)+Y*C(3) > 

TEMP«TEMP/<1.  D0+Y*<D<1)+Y*<D<2>+Y*D<3)>  >> 
0AU8I V-SNOL ( TEMP-Y ) 

IF ( Z.  OT.  0.  9D0)  OAUSIV— OAUSIV 

RETURN 

END 


A.  1.13  Weibull  Background  Generator  (WEIBUL  and  SETUPW) 

Calling  Statement ;  Call  WEIBUL  (XR,  XI,  N,  IREP) 

Parameters  Read  from  Input  File; 

ISEEDC  ■  starting  seed  for  pseudo-random  number  generator 
C  *  exponent  parameter  of  Weibull  clutter  power 
distribution 

Argument  List  Definition: 

Input  Variables:  N  ■  number  of  complex  samples 

generated  by  each  call 

IREP  -  representation  code 

Returned  Variables:  XR  and  XI  are  defined  according 

to  IREP,  as  discussed  below. 

Other  Modules  Called:  SETUPW,  RANAMP,  RANPHS 
Description : 

For  IREP  *  -1,  WEIBUL  generates  N  uncorrelated  complex  samples 
of  Weibull  video,  with  unit  mean  power  and  uniformly  distributed  phase. 

That  is 

XR  -  /T  sin  $ 

XI  -  /T  cos  $ 

where  Y  has  Weibull  distribution  with  probability  density 

C— 1  c 

p(Y)dy  -  —  (^-)  exp  {-(—)  }  dy  ,  Y  ^  0 
o  o  o 


C  *  exponent  parameter 

Y  -  1.0/r(l  +  1/C) 
o 


where 


and  where 


<p  “  Can  *  (XR,XI)  is  uniform  on  (-ir.ir). 

The  method  of  generation  is  identical  to  that  used  in  the  log- 
normal  generator  LGNORM,  except  that  7  has  a  Weibull  distribution.  The 
inverse  probability  table  is  generated  with  a  call  to  SETtIPW.  . 


SUBROUTINE  WEXBUL(XR,  XI*  N* IREP) 

C  WE I BULL  CLUTTER  GENERATOR 

C 

C  I REP  RESULT 

C  0:  RETURNS  NATURAL  LOGARITHM  OF  AMPLITUDES 

C  IF  IREP  WERE  1  (SEE  BELOW). 

C 

C  1:  RETURNS  N  UNCORRELATED  WEIBUL  AMPLITUDES 

C  IN  XR.  WITH  UNIT  MEAN  POWER.  AND  EXPONENT 

C  PARAMETER  C  (POWER). 

C 

C  2:  RETURNS  N  UNCORRELATED  WEIBUL  POWERS  IN  XR. 

C 

C  -1:  RETURNS  N  UNCORRELATED  WE I BULL  CLUTTER  SAMPLES 
C  WITH  RANDOM  PHASE.  UNIT  MEAN  POWER.  AND 

C  EXPONENT  PARAMETER  C.  RETURNS  NSAMPLES  IN  XR. 

C  AND  N  QUADRATURE  SAMPLES  IN  XI.  USES  3-LEVEL 

C  TABLE  LOOK-UP  FOR  AMPLITUDES.  UNIT  CIRCLE  ALGORITHM 

C  FOR  SIN  AND  COS. 

C 

C  -2:  RETURNS  N  SIN  Si  COS  PAIRS  ONLY.  IN  XR  Si  XI 
C 

C  -3:  RETURNS  N  SCN(SIN)  Si  SGN(COS)  PAIRS  ONLY 
C 

•INSERT  JIMB  AN I T A>SUR ADS>  IOLUS 
•INSERT  JIMB  AN I T A>SUR ADS>  FLAGS 
INTEGER*4  ISEEDC.  I SEED 
DIMENSION  XR(1>, XI(1). PIV(3S4> 

DATA  PIV.  LPI V.  NPIV.  ISEED/384*1.  O.  3.  7.  323889239/ 

C 

C  ON  FIRST  CALL  READ  PARAMETERS  AND  INITIALIZE  PIV  TABLE 
C 

IF(.  NOT.  ( IFLAGW.  EQ.  0)  )  GO  TO  29 
C 

WRITE ( LUOUT .  100) 

100  F0RMAT(////1X. 'WE I BULL  CLUTTER  GENERATOR') 

READ(LUIN.  »)  ISEFDC 
IF< ISEEDC.  EQ.  O)  ISEEDC -I SEED 
WRITE (LUOUT.  200)  ISEEDC 
200  FORMAT (//AX.  'PRN  GENERATOR  SEED  -  ',112) 

READ(LUIN,  *)  C 
WRITE (LUOUT,  300)  C 

300  FORMAT ( //6X»  'EXPONENT  PARAMETER  -  ',F7.  3) 

WRITE (LUOUT.  400)  IREP 

400  FORMAT (//AX.  'REPRESENTATION  CODE  »  M3) 


CALL  SETUPW(PIV» NPIV,  LPIV,  IREP, C> 

IFLAGW" 1 

CONTINUE 

IF(  IREP.  LT.  -1 )  GO  TO  28 

GENERATE  RANDOM  AMPLITUDES  (OR  POWERS) 

CALL  RANAMP(PIV»  NPIV.  LPIV, ISEEDC,  XR,  N) 

IF  IREPCO  GENERATE  AND  APPLY  RANDOM  PHASE 
CONTINUE 

IF( IREP.  OE.  0)  GO  TO  30 

CALL  RANPHS( XR, XI,  N.  ISEEDC.  IREP) 

RETURN 

END 


■  ■  ■*)£■' Jk. 


f 


Of»f> 


SUBROUTINE  SETUPUtPIV.  NPIV,  LPIV,  ZREP.C) 

SET  UP  INVERSE  PROBABILITY  TABLE  FOR  WE I BULL  GENERATOR. 

DOUBLE  PRECISION  PROB 
♦INSERT  JIMB  AN I T A>SUR ADS>  IOLUS 
DIMENSION  PIV(1> 

C 

CINV-1.  O/C 
CONST-1.  O+CINV 
CONST-1.  /GAMMA <  CONST  > 

IF<IREP.  EG.  2)  GO  TO  9 
CINV-CINV*.  9 
CONST-SORT < CONST > 

9  CONTINUE 

NN-2**NPIV 
JLP— 0 

DO  20  LP— 1,  LPIV 
DO  10  J-l.NN 
JLP-JLP+1 

PROB— 1.  DO+DBLE( FLOAT <J-NN>-0.  9) /( DBLE < FLOAT ( NN ) )**LP) 
PIV( JLP )  — (SNGL(-DLOG( 1.  DO-PROB ) ) **C I NV > *CONST 
IFCIREP. EQ.  O)  PIV( JLP )-ALOG<PIV< JLP > ) 

CONTINUE 
CONTINUE 


10 

20 

C 


CALL  WRTPIV<PIV,  NPIV,  LPIV,  LUOUT) 

RETURN 

END 
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A. 1.14  Gamma  Function  (GAMMA) 


Calling  Statement:  Y  *  GAMMA(X) 

Description : 

GAMMA  computes  the  gaama  function,  using  the  fact  that 

T(x  +  1)  -  xr(x) 

and  a  polynomial  approximation  for  1.0  <  x  <  2.0  from  equation  6.1.35  of  [7]. 
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c 

c 

c 

c 


10 

c 

c 

c 


20 

c 

c 

c 


30 

C 

C 

C 


FUNCTION  GAMMA(Y> 

GAMMA  FUNCTION  POLYNOMIAL  APPROXIMATION 
REF.:  ABRAMOWITZ  *  STEGUN  P.  257.  EQ.  6.  1. 35 

DIMENSION  A(5> 

DATA  A/-.  5748646.  .  9512363.  -.  6998588.  .  4245549.  1010678/ 

X-Y 

GAMMA" 1 . 0 
CONTINUE 

IF<X.  LE.  2.  0)  GO  TO  20 

X>2.  0 

X-X-l.  0 
GAMMA-GAMMA* X 
GO  TO  10 
CONTINUE 

IFCX.  GE.  1.  O)  GO  TO  30 
X<1.  O 

GAMMA-GAMMA/ X 
X-X+l.  O 
GO  TO  20 
CONTINUE 

1.  0C-XO2.  O 

X-X-l.  O 

GAMMA— GAMMA* ( 1 .  0+< A<  1  >  +  < A<2>  +  < A(3)  +  ( A<4>+A(5)*X >*X )*X )*X>*X ) 

RETURN 

END 


A.  1.15  Correlated  Background  Generator  (CORBKG) 


Calling  Statement;  Call  COBBKG  (XR,  XI,  XMOD,  N,  K,  IREP,  ICMBIN) 

Parameters  Read  from  Input  File: 

ITYPE  -  type  of  fluctuations  for  average  cross-section 
I TYPE  »  2  “  log-normal 
ITYPE  -  3  -  Weibull 

Argument  List  Definition; 

XR,  XI,  XMOD  -  are  components  of  background  samples  returned  by 
CORBKG ,  defined  according  to  IREP  and  ICMBIN, 
as  discussed  below. 

N  *  number  of  range  cells 

K  -  number  of  pulses  per  range  cell 

IREP  *  representation  code 

ICMBIN  =  combination  code  to  determine  whether  components 
will  be  combined  (see  discussion  below) . 

Other  Modules  Called:  NOISE,  LGNORM,  WEIBUL 


Discussion : 

CORBKG  generates  log-normal  or  Weibul  background  samples  which  are 
Rayleigh  moduled  pulse- to-pulse.  Samples  are  generated  on  each  of  N  range 
bings,  and  uncorrelated  modulations  are  generated  for  each  bin  for  K  pulses. 
The  sample  are  generated  with  a  call  to  either  LGNORM  or  WEIBUL.  The  Rayleigh 
fluctuations  are  generated  by  NOISE. 

Sample  representations  are  determined  by  IREP,  just  as  they  are  for 
the  NOISE,  LGNORM,  and  WEIBUL  generators.  Each  sample  is  the  product  of  two 
components,  the  log-normal  or  Weibull  average  cross-section,  and  the  Rayleigh 
fluctuation  about  that  average  cross-section.  If  ICMBIN  -  0,  the  components 
are  returned  separately.  If  ICMBIN  ■  1,  they  are  combined  prior  to  return. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  COR8KG(XR.  XI.  XMOD.  N. K,  IREP.  ICMBIN) 

GENERATES  LOO-NORMAL  OR  WEIBULL  BACKGROUND  SAMPLES  WHICH 
ARE  RAYLEIGH  MODULATED  PULSE-TO-PULSE.  MEAN  POWER  OF 
RESULT  IS  UNITY.  SEE  W.  RIVERS.  “DISTRIBUTION  FUNCTION 
OF  A  RAYLEIGH  PROCESS  MODULATED  BY  A  LOG-NORMAL  PROCESS'S 
TSC-W25— 99.  N  MODULATIONS  OF  MEAN  POWERS  ARE  RETURNED  IN  XMOD. 
N*K  RAYLEIGH  FLUCTUATIONS  ARE  RETURNED  IN  XR  AND  XI. 


IF  ICMBIN-0 

IREP  XR 

2  RAYLEIGH 

POWER 

1  RAYLEIGH 

AMPLITUDE 

0  LOO (RAYLEIGH 
AMPLITUDE) 

-1  RAYLEIGH  I 

-2  B.  P.  LIM.  I 

-3  HARD-LIM.  I 


RESULTS  ARE  AS  FOLLOWS: 


XI 


RAYLEIGH  G 
B.  P.  LIM.  Q 
HARD-LIM. G 


XMOD 

LOG-NORM  OR  WEIBULL  POWER 
LOG-NORM  OR  WEIBULL  AMPLITUDE 
LOG ( LOG-NORM  OR  WEIBULL  AMPLITUDE) 
LOG-NORM  OR  WEIBULL  AMP. 


IF  ICMBIN-1.  MODULATIONS  ARE  APPLIED  TO  XR  AND  XI. 


9 INSERT  JIMB  ANITA>SURADS>IOLUS 
9 INSERT  JIMB  AN I T A>SURADS>FLAOS 

DIMENSION  XR(1). X I ( 1 ) . XMOD ( 1 ) 
DATA  I COR/2/ 


C 

C 

C 


100 


200 

C 

29 


C 

C 

C 


90 

900 


ON  FIRST  CALL  INITIALIZE 

IFC.NOT.  (IFLAG(ICOR).  EG.  0))  GO  TO  29 
WRITE(LUOUT.  100) 

FORMAT (////IX.  'CORRELATED  BACKGROUND  GENERATOR 
READ(LUIN.  *)  I TYPE 
WRITE (LUOUT.  200)  I TYPE 

FORMAT ( //6X. 'SPATIAL  MODULATION  TYPE  - 
I FLAG ( ICQR)*1 


13) 


CONTINUE 

NBIN“N*K 

IREPP-IABS(IREP) 

CALL  NOISE(XR. XI.NBIN.  IREP) 

IF( IREP.  LT.  -1 )  GO  TO  900 

IF ( I TYPE.  EG.  2)  CALL  LGNORM( XMOD. XMOD. N, IREPP) 
IF(  ITYPE.  EG.  3)  CALL  WEIBUL (XMOD.  XMOD.  N,  IREPP) 
IF( ICMBIN.  EG.  0)  GO  TO  900 

COMBINE  COMPONENTS  APPROPRIATELY 

IF(IREP.  LT. -1)  GO  TO  900 
JJ-1 

DO  90  KK-l.NBIN 
XMODJJ-XMOD(JJ) 

IF (IREP.  NE.  0)  XR ( KK ) “XR ( KK ) *XMODJJ 
IF ( IREP.  EG.  -1)  X I ( KK ) "X I ( KK > *XMODJ J 
IF ( IREP.  EG.  0)  XR ( KK ) *XR ( KK ) +XMOD J J 
JJ-JJ+1 

IF(JJ.  OT.  N>  JJd 
CONTINUE 
RETURN 
END 
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A.  1.16  Radar  IF  Generator  (IFGEN) 


Calling  Statement:  Call  IFGEN  (YR,  YI,  XR,XI,  NXY,  HXY,  W,  NW) 

Argument  List  Definitions: 

Input  Variables:  XR,XI  -  matrix  of  input  signals 

MXY ,  NXY  »  row  and  column  dimensions 
of  XR,  XI,  YR,  and  YI 

W  “  vector  of  waveform  sample 

NW  »  length  of  W 

Output  Variables:  YR,  YI  ■  matrix  of  output  signals 

Description: 

IFGEN  computes  the  convolution  of  the  rows  of  the  complex  signal 
X  with  the  real  function  W.  It  is  assumed  that  all  samples  of  W  are  either 
+1  or  -1.  There  is  no  wrap-around  on  the  ends.  Thus,  only  the  first 
NXY  -  NW  +  1  columns  of  each  row  of  Y  are  computed.  The  output  YR  is 
equivalent  to: 


NW 

YR(IY+( JY-1) *NXY)  -  Z  W(IW)*IR(IY+NW-IW+(JY-1)*NXY) 

IW-1 


for 

IY  -  1,2,..., NXY  -  NW  +1 
JY  -  1,2,..., MXY 

YI  is  similarly  defined  in  terms  of  XI. 

This  subroutine  is  used  to  generate  the  radar  IF  signal  (actually 
bipolar  video)  by  convolving  the  transmit  waveform  W  with  the  distributed 
clutter  samples  XR,  XI. 
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SUBROUTINE  IFGEN  (  YR,  YI,  XR,  XI,  NXY,  MXY,  W,  NW> 

COMPUTES  CONVOLUTION  OF  COMPLEX  SIGNAL  X  WITH  REAL 
FUNCTION  W.  WHERE  EACH  SAMPLE  OF  W  IS  +/-  1.0. 
RESULT  IS  RETURNED  IN  Y. 

SHAPE  OF  X  IS  MXY  ROWS  BY  NXY  COLUMNS*  STORED 
AS  A  ONE  DIMENSIQNSL  ARRAY.  ROW  BY  ROW. 

SHAPE  OF  Y  IS  THE  SAME  AS  X. 

THE  CONVOLUTION  IS  WITHIN  EACH  ROW  OF  X.  FOR  THE 
FIRST  ROW.  THE  RESULTS  ARE  EQUIVALENT  TO: 


Y(IY)«SUM(  W< IW)*X < IY+NW-IW)  ) 

IW-1 

FOR  IY-1.NY 

WHERE  NY-NX-NW+1 

THUS,  ONLY  THE  FIRST  NY  COLUMNS  OF  EACH  ROW  OF 
Y  CONTAIN  RESULTS. 

CAUTION:  X&Y  MUST  BE  DISTINCT  ARRAYS  IN  THE  CALLING 
PROGRAM 

DIMENSION  YR < 1 ) .  YK1J,  XR(1>,  XI(1).  W(l) 

NY-NXY-NW+1 

IK-1 

DO  20  K-l.MXY 

DO  10  I-l.NXY 

YR ( IK)— O.  O 
YI(IK)-0.  O 
IK-IK+1 
CONTINUE 
CONTINUE 

DO  80  IW-l.NW 

KOFFST— NXY 

IF  (.NOT.  (W(IW).  GT.  0.  0))  GO  TO  90 
DO  40  K-l.MXY 

KOFFST -KOFFST +NX  Y 
I X -NW+ 1 - I W+KOFFST 
DO  30  I-l.NY 

I Y— KOFFST ♦ I 

YR< IY)— YR ( IY)+XR ( IX ) 

YI(IY>«YI(IY)+XI(IX> 

IX-IX+1 

CONTINUE 

CONTINUE 

CONTINUE 

IF  (.NOT.  (W(  IW) .  LT.  0.  0)  )  GO  TO  80 
DO  70  K-l.MXY 

KOFFST-KOFFST+NXY 
IX— NW+1- I W+KOFFST 
DO  60  I-l.NY 

IY— KOFFST+I 

YK(IY)-YR(IY)-XR(IX) 

YKIY)-YKIY)-XI(IX) 

IX-IX+1 

CONTINUE 

CONTINUE 

CONTINUE 

RETURN 

END 
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A.  1.17  Target  Generator  (TARGET) 


Calling  Statement:  Call  TARGET  (X,  N,  K) 

Parameters  Read  from  Input  File; 

XSEEDT  -  starting  seed  for  pseudo-random  number  generator 

ISWRL  -  Swerling  Case  number  (0,1, 2,3,4) 

Argument  List  Definitions: 

Input  Variables:  N  -  number  of  independent  target 

samples  per  pulse 

K  *  number  of  pulses 

Returned  Variables:  X  ■  vector  of  NxK  target  samples 

Other  Modules  Called:  CHISQ,  WRTPIV,  RANAMP 

Description: 

TARGET  generates  N  x  K  target  amplitude  samples  of  unit  mean 
power.  N  is  the  number  of  target  samples  per  pulse,  and  K  is  the  number  of 
pulses.  Target  distributions  correspond  to  Swerling  cases  0-4.  For  case 
0,  X  contains  N  x  K  ones.  For  cases  2  and  4,  X  contains  N  x  K  uncorrelated 
samples.  For  cases  1  and  3,  X  contains  N  uncorrelated  samples,  repeated 
K  times  in  blocks  of  N  samples. 

The  target  amplitudes  are  given  by 
X  -  /T 


where  Y  is  a  unit  mean  random  variable  whose  probability  density  function 
is  given  by 


py(y) 


S(Y-1.0)dy,  where  5 (Y) 
function, 
case  0, 


is  the  delta 
for  Swerling 


exp(-Y)dy,  Swerling  case  1  & 
4Yexp(-2Y)dy,  Swerling  cases 


2,  y  >  0 
3  &  4,  y  _>  0 


The  method  of  generation 
and  table  look-up. 


is  by  inverse  probability  table  construction 
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A.  1.17  Target  Generator  (TARGET) 


Calling  Statement:  Call  TARGET  (X,  N,  K) 


Parameters  Read  from  Input  File; 


ISEEDT  ■  starting  seed  for  pseudo-random  number  generator 
ISWRL  -  Swerling  Case  number  (0,1, 2, 3, 4) 


Argument  List  Definitions: 


Input  Variables: 


Returned  Variables: 


N  -  number  of  independent  target 
samples  per  pulse 

X  »  number  of  pulses 

X  -  vector  of  N  x  X  target  samples 


Other  Modules  Called:  CHISQ,  WRTPIV,  RAN AMP 
Description: 

TARGET  generates  N  x  K  target  amplitude  samples  of  unit  mean 
power.  N  is  the  number  of  target  samples  per  pulse,  and  K  is  the  number  of 
pulses.  Target  distributions  correspond  to  Swerling  cases  0-4.  For  case 
0,  X  contains  N  x  K  ones.  For  cases  2  and  4,  X  contains  N  x  X  uncorrelated 
samples.  For  cases  1  and  3,  X  contains  N  uncorrelated  samples,  repeated 
K  times  in  blocks  of  N  samples. 

The  target  amplitudes  are  given  by 

X  -  ST 

where  Y  is  a  unit  mean  random  variable  whose  probability  density  function 
is  given  by 


Py(Y) 


6(Y-1.0)dy,  where  5(Y)  is  the  delta 
function,  for  Swerling 
case  0, 

exp(-Y)dy,  Swerling  case  1  &  2,  y  ^  0 
4Yexp(-2Y)dy,  Swerling  cases  3  &  4,  y  ^  0 


The  method  of  generation  is  by  inverse  probability  table  construction 
and  table  look-up.’ 
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SUBROUTINE  TARGET (X.  N.  K) 

C 

C  GENERATES  N*K  TARGET  AMPLITUDE  SAMPLE3  OF  UNIT  MEAN 

C  POWER.  N  IS  NUMBER  OF  RANGE  BINS.  K  IS  NUMBER  OF 

C  PULSES.  GENERATES  SAMPLES  FOR  SUER LI NO  CASES  0-4. 

C  RETURNS  SAMPLES  IN  X.  FOR  CASE  0.  X  CONTAINS 

C  N*K  ONES.  FOR  CA8ES  2&4.  X  CANTAINS  N«K  UNCORRELATED 
C  SAMPLES.  FOR  CASES  1*3.  X  CONTAINS  N  UNCORRELATED 

C  SAMPLES.  REPEATED  K  TIMES  IN  BLOCKS  OF  N  SAMPLES. 

C 

INTEGER*4  ISEEDT.  I SEED 
DOUBLE  PRECISION  PROB. DELPRB.  PO 
•INSERT  IOLUS 
•INSERT  FLA08 

DIMENSION  X(1).PIV(384) 

DATA  PIV.  LPIV.  NPIV.  XSE£D/384«1.  0.  3.  7.  924922231/ 

C 

C  ON  FIRST  CALL  READ  PARAMETERS  AND  INITIALIZE  PIV  TABLE. 

C 

IF  (.NOT.  (IFLAOT.  EG.  0) )  GO  TO  29 
URITE(LUOUT.  100) 

100  FORMAT <//// IX.  'TARGET  AMPLITUDE  GENERATOR') 

READ(LUIN< *)  ISEEDT 
IF < ISEEDT.  EG.  0)  ISEEDT* I SEED 
WRITE(LUOUT.  200)  ISEEDT 
200  FORMAT (//6X» 'PRN  GENERATOR  SEED  *',I12> 

READ (LU IN,  *)  ISWRL 
WRITE(LUOUT,  300)  ISWRL 

300  FORMAT! //6X.  'SWEKLING  CASE  '.13.  '  TARGET  MODEL') 

C 

M*( ISWRL+1 )/2 
PO-O.  O 
NN-2**NPIV 
JLP*0 

DO  20  LP-l.LPIV 

DELPRB*.  29D0/ (DBLE( FLOAT <NN> )**LP) 

DO  10  J*1.NN 
JLP-JLP+1 

PR0B*1 .  DO+DBLEC FLOAT (J-NN)-O.  9)/(DBLE(FLQAT(NN> >**LP> 
PIVJLP*CHISQ(PROB.  DELPRB.  M.  PO) 

IF(M.  NE.  0)  P I V JL:  «P I V  JLP /FLOAT ( M ) 

P IV  <  JLP ) -SORT ( PIVJLP ) 

10  CONTINUE 

20  CONTINUE 

IFLAOT- 1 

CALL  WRTP IV  <  P IV.  NH I V,  LP IV.  LUOUT ) 

C 

29  CONTINUE 
C 

C  GENERATE  RANDOM  AMPLITUDES 

C 

NSAMP-N 

IF(MOD< ISWRL. 2).  EG.  0)  NSAMP*N*K 

CALL  RANAMP (PIV,  NPIV,  LP IV.  ISEEDT. X,  NSAMP) 

IF ( MOD (  ISWRL,  2).  EG.  0)  GO  TO  90 
IF(K.  EQ.  1)  GO  TO  90 
DO  40  1*1. N 
XX*X(I) 

IJ*I 
KM1-K-1 
DO  30  J-l,  KM1 
IJ-IJ+N 
X(IJ)-XX 

30  CONTINUE 
40  CONTINUE 

C 

90  RETURN 
END 
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A. 1.18  Inverse  Cumulative  Chi-Squared  Function  (CHISQ) 


Calling  Statement:  T  ■  CHISQ (F,  DP,  M,  TOLD) 

Argument  List  Definition: 

Input  Variablea :  P  -  probability  level  (double  precision) 

D?  -  probability  tolerance  desired 
(double  precision) 

M  “  nuaber  of  duo -degrees  of  freedom 
(i.e.  nuaber  of  degrees  of 
freedom  *  2  xX) 

TOLD  ■  previous  value  returned  by  CHISQ 
(dduble  precision;  initially 
zero) 

Returned  Variables:  CHISQ  -  aaplltude  for  probability 

level  P  (single  precision) 

TOLD  ■  double  precision  replica 
of  CHISQ. 


Other  Modules  Called:  GAUSIV 
Description: 

CHISQ  returns  the  Inverse  cumulative  probability  distribution 
function  for  a  Chi-squared  random  variable  vith  M  duo-degrees  of  freedom* 
and  mast  value  M.  For  a  given  probability  level  P,  the  P-level  amplitude 
is  obtained  by  inverting  the  cumulative  distribution  function  given  by. 

M-l  _k 

P(T)  -  1  -  aacp(-T)  Z  rr  ,  M  >  1 
k-0 

For  H«0,  a  value  of  CHISQ  -  1.0  is  always  returned. 

The  inversion  of  P(T)  is  obtained  in  the  following  manner. 

For  M  •  1, 

T(P>  -  -  ln(l  -  P). 
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for  1  <  M  <  50,  1(f)  la  found  by  unking  an  Initial  guass  Yp*  and  counting 
fClp).  If  PCIp)  la  within  DP  of  P,  tha  Inversion  la  dona,  and 

CH1SQ  -  Tp 

la  raturaad.  If  not,  Chan  a  naw  guaaa  la  obtainad  aa 

P  -  P(Y_) 

V  “  YP  +  3P (Tp) 

3T 

and  tha  proeaaa  rapaatad  until  P(Yp  )  la  eloaa  enough  to  P.  Tha  Initial 
guaaa  for  Yp  la  found  by  althar  of  two  tachnlquaa.  On  tha  first  call  to 
CHISQ  (signified  by  TOLD  ■  0) ,  Yp  is  obtainad  by  sealing  tha  M  -  1  root 
using  a  Gaussian  approximation: 

>  lIL  .  “la(l  -  P)-l 

& 

On  subsaquant  calls,  tha  previous  root  la  used: 

Tp  -  TOLD 

For  M  >50,  tha  Chi-Squara  distribution  is  approximated  as  Gaussian, 
and  tha  inversion  is  obtained  by  calling  GADS IV. 
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ooo  non  non  n  n  nnnnnnnn 


REAL  FUNCTION  CHISQ  (PROS.  DELPRB*  M* XOLO) 

INVERSE  PDF  FOR  CHI-SQUARE 

DEGREE  2«M.  MEAN-M.  RANDOM  VARIABLES. 

RETURNS  CHISQ-1.0  FOR  M-O. 

USES  OAUSSIAN  APPROX.  FOR  M>90. 

DOUBLE  PRECISION  PROBABILITY  VERSION 

DOUBLE  PRECISION  PROB.  DELPRB,  XOLD* X.  XM*  PX.  EXPO 
DOUBLE  PRECISION  XTOK.  SUM 

CHISQ-1.  0 

IF(M.  EQ.  0)  RETURN 

IF  (M  .  OE.  91)  00  TO  200 

X  -  -DL0G(1.D0  -  PROB) 

FOR  M-l.  THIS  IS  AN  EXACT  SOLUTION 

IF  (M  .EQ.  1)  00  TO  100 

IF  (.NOT.  (XOLD  .EQ.  0.0))  00  TO  9 

SCALE  M— 1  ROOT  ON  OAUSSIAN  APPROXIMATION 

XM  -  DBLE ( FLOAT ( M ) ) 

X  -  XM  -  DSQRT ( XM ) * < 1 .  DO  -  X) 

00  TO  10 

START  AT  PREVIOUS  ROOT  LOCATION 

X  -  XOLD 
CONTINUE 
XTOK  -  1.  DO 
SUM  -  0.  0 
DO  30  I  ■  l.M 
K  -  I  -  1 

IF  (K  .EQ.  0)  00  TO  20 

XTOK  -  XTOK  *  X/DBLE ( FLOAT (K)) 
8UM  -  SUM  +  XTOK 
CONTINUE 
EXPO  -  DEXP(-X) 

PX  -  1.  DO  -  EXPO  #  SUM 

IF  (DABSCPRQB  -  PX)  .  LE.  DELPRB)  00  TO  100 
X  -  X  ♦  (PROB  -  PX)/(EXPO  *  XTOK) 

00  TO  10 
CONTINUE 
CHI SO  -  SNOL(X) 

XOLD  -  X 
RETURN 

CONTINUE 

IF  (PROB  .OT.  .  9D0)  00  TO  210 
Y  -  0AUS1 V ( PROB  > 

CHI8Q  -  FLOAT(M)  +  Y*SQRT (FLOAT ( M) ) 

RETURN 

CONTINUE 

Y  •  OAUSIVd.  DO  -  PROB) 

CHISQ  -  FLOAT (M)  -  Y*SGRT  <  FLOAT  CM)) 

IF  (CHISQ  .LT.  0.0)  CHISQ  -0.0 

RETURN 


A.  1.19  Threshold  Crossing  Accuaulatlon  (PDDFX) 
Calling  Stateaent:  Call  PDDFXQO 


Argument  List  Definition: 

M  ■  number  of  samples  in  waveform  WR  (in  labelled  common  ARRAYS) 

Description: 

This  is  the  module  which  actually  implements  the  receiver  signal 
processing.  The  IP  background  signals  are  presumed  to  be  in  array  VUR  and 
VUI.  The  target  amplitudes  are  in  the  array  TGT.  The  waveform  used  in 
generating  VUR  and  VUI  is  in  the  array  WR.  It  is  assumed  that  WR  is  either 
+1  or  -1. 

PDDFX  generates  a  target  waveform  and  adds  it  to  the  background 

IF  signal  at  the  appropriate  signal-to-background  ratio.  The  in-phase  and 

quadrature  sum  signal  are  than  hard-limited  to  +/-1  and  pulse  compressed 

using  the  waveform  WR.  A  test  statistic  (TEST)  if  formed  by  summing  the 
2  2 

pulse- compressed  I  +  Q  over  NPDI  pulses.  The  test  statistic  is  compared 
to  the  array  of  thresholds, and  the  threshold  crossing  counts  accumulated  in 
the  matrix  NIC.  The  test  statistic  and  threshold  comparison  is  repeated  on 
each  range  bln.  The  process  is  then  repeated  for  all  other  signal-to- 
background  ratios  SBR.  On  return,  the  matrix  NTC  contains  the  cumulative 
count  of  threshold  crossings  as  a  function  of  threshold  and  signal-to- 
background  ratio. 
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SUBROUTINE  PDDFX(M) 

C 

C  COMPUTE  DETECTION  PROBABILITY  FOR 

C  DICKE-FIX  CFAR  (HARD-LIMITED  I  *  Q). 

C 

C  ADD  TARGET  PLUS  BACKGROUND.  PERFORM  HARDLIMITINO 
C  AND  PULSE  COMPRESSION.  SUM  PULSES  INCOHERENTLY. 

C  AND  TEST  AGAINST  THRESHOLDS. 

C 

*  INSERT  J.IMB  JODI>SURADS>CNTROL 
•INSERT  JIMB  J0DI>SURAD8>I0LUS 
•INSERT  JIMB  JOD I >SURAD8> ARRAYS 
MH1T2— (M— 1 )*2 
NRM2-NRBIN+MM1T2 
DO  20  J— 1.  N8BR 
SB-SBR(J) 

I8BR-  <  J— 1 )  -*NTHR 
DO  19  K-l.NRBIN 
TEST-0.  0 
KL-K-NRM2 
KLT-K-NKBIN 
DO  9  L-l.NPDI 
KL-KL+NRM2 
KLT-KLT+NRB IN 
KLM-KL 
TESTR-O.  0 
TEST I -0.  0 
TARG='SB*TGT  <  KLT  > 

DO  4  MM— 1 »  M 

IF  (.NOT.  (MR  (MM).  GT.  0.  0)  )  00  TO  31 

TESTR— TE8TR+SIGN( 1 .  O.  ( VUR ( KLM ) +TARG ) ) 
TEST I -TEST I +S I ON < 1.  0.  VUI (KLM) ) 

GO  TO  41 

31  CONTINUE 

TESTR— TESTR-S I ON ( 1 .  0.  (VUR(KLM)— TARG) ) 
TEST I -TEST I -8 1 GN(1.  0.  VUI (KLM)) 

41  CONTINUE 

KLM— KLM+1 

4  CONTINUE 

TESTI  -TE8TR*TESTR+TEST I *TEST I 
TEST— TEST +TESTT 

9  CONTINUE 
INDX-ISBR 

DO  10  IT-l.NTHR 
INDX-INDX+1 

IF(Tfc8T.  OT.  THR(IT) )  NTC ( INDX )-NTC < INDX >+l 

10  CONTINUE 

19  CONTINUE 


20  CONTINUE 

RETURN 
END 
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A.  1.20  Thru  hold  Crossing  S  usury  (SUMARY) 


ii  Call  SUMAKY  (NIC,  SBR,  THR,  HSBR,  HTHR,  HTRIAL) 


Input  Variables: 


HTC  »  array  of  thruhold  crossing 
counts  (double  precision) 

SBR  -  array  of  si goal- to-bsckgro  und 
ratios 

TER  •  array  of  thresholds 
HSBR  -  nuaber  of  elements  in  SBR 
NTER  -  nuaber  of  elements  in  TER 
HTRIAL  •  nuaber  of  trials  accumulated. 


Returned  Variables : 


SUMAKY  converts  the  threshold  crossing  counts  stored  in  HTC  to 
threshold  crossing  probabilities ,  and  vrltu  a  sum ary  table  to  logical 
unit  LUOIJT  (passed  via  labelled  coamon  IOUJS) .  The  array  HTC  is  a  single 
dimension,  double  length  integer  array  with  crossing  counts  stored  in  the  first 
HSBR  x  NZHR  locations.  The  first  HTHR  locations  contain  counts  for  the 
first  value  of  SBR.  The  next  HTHR  locations  are  for  the  second  value  of 
SBR,  and  so  forth. 


c 

c 

c 

c 


SUBROUTINE  SUMARY (NTC.  SBR,  THR.  NSBR.  NTHR,  NTRIAL) 

SUMMARIZES  THRESHOLD  CROSSING  COUNTS 
RECORDED  IN  NTC. 


DOUBLE  PRECISION  XTRIAL.  PROS 
INTEGER *4  NTRIAL.  NTC 
•INSERT  IOLUS 

DIMENSION  NTC<1).PR0B(20).SBR(1).  THR<1) 
DATA  NCOL/4/ 


C 

100 

200 


300 

400 


10 

900 

20 


30 

C 


WRITEtLUOUT.  100) 

F0RMAT(////1X.  'SUMMARY  OF  THRESHOLD  CROSSING  PROBABILITIES(PCT.  )') 
WRITE(LUOUT>  200)  NTRIAL 


FORMAT < //AX.  'NO.  OF  TRIALS  - 


17) 


10611.  3/) 


XTRIAL-100.  DO/DBLE< FLOAT (NTRIAL) ) 
NLOOP-NTHR/NCOL 
nnloop ■NCOL*NLOOP 
IF?NNLOOPMIt.  NTHR >  NLODP-NLOOP+1 
Jl-1 

J2*J1 +NC0L-1 
DO  30  I-l.NLOOP 
J2«MINO( J2,  NTHR ) 

WRITE(LUOUT«  300) (THR(U). J-Jl,  J2) 

FORMAT ( ///AX.  'THRESHOLD  CONSTANT 
URITE(LUOUT.  400) 

FORMAT ( /6X»  'S/B  RATIO(DB)  ') 

DO  20  K-1.N8BR 
DBSBR-SBR(K) 

I SBR* ( K— 1 ) *  NTHR 
DO  10  J-J1.J2 
INDX-ISBR+J 

PROS  <  J ) -DBLE ( FLOAT ( NTC (INDX)))*XTRIAL 
CONTINUE 

WRITE<LUOUT» 900)  DBSBR.  CPROB< J).  J-Jl,  J2) 
FORMAT ( /AX.  FI 0.  1.  11X.  9F11.  9) 

CONTINUE 

J1-J1+NC0L 

J2-J2+NC0L 

CONTINUE 

RETURN 

END 
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A.  1.21  Edit,  Fit,  and  Plot  Data  (PLOTPD)' 


Calling  Statement;  Call  PLOTPD  (NTC,  SBR,  NSBR,  NTHR,  NTRIAL) 

Parameters  Read  from  Input: 

IFNAM  -  plot  file  name 

PFAPCT  ■  vector  of  threshold  crossing  probabilities  for  no  signal. 

Argument  List  Definition: 

NTC  •  matrix  of  threshold  crossing  counts 

SBR  ■  vector  of  signal- to-background  ratios  In  dB. 

NSBR  -  length  of  SBR 

NTHR  ■  number  of  thresholds  summarized  in  NTC 

NTRIAL  -  number  of  trials 

Other  Modules  Called:  PDPLOT,  SEARCH 

•  Description : 

PLOTPD  First  opens  a  plotfile  for  writing  plottable  curves 
of  PQ  vs.  SBR.  For  each  threshold  level,  PLOTPD  then  extracts  the  appropriate 
threshold  crossing  counts  from  NTC,  edits  them,  and  converts  to  percent 
probability  of  threshold  crossing.  These  probabilities  are  then  distorted 
using  the  value  of  PFAPCT.  The  distorted  values  are  passed  to  PDPLOT  for 
curve  fitting  and  plot  file  writing. 

The  data  editing  rejects  crossing  counts  of  fewer  than  10  events, 
or  greater  than  10  less  than  the  total  number  of  trials.  In  addition, 
it  rejects  crossing  probabilities  less  than  twice  the  false  alarm  probability. 
Finally,  it  rejects  probabilities  which  are  less  than  a  probability  obtained 
at  a  lower  signal-to-background  ratio,  i.e.  it  assures  that  the  probabilities 
passed  to  PDPLOT  will  be  monotomlcally  increasing  with  signal-to-background. 

The  probability  distortions  are  to  enable  curve  fit  prameters  to 
be  obtained  using  the  approximate  procedure  discussed  in  ref  [2J .  The 
distortion  is: 


where  all  probabilities  are  in  percent. 


-  tjl 
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SUBROUTINE  PLOTPD  (NTC,  SBR, NSBR, NTHR,  NTRIAL) 

CONVERTS  THRESHOLD  CROSSING  COUNTS'  IN  NTC  TO 
DETECTION  PROBABILITIES  IN  PERCENT.  ELIHINATES 
STATISTICALLY  INSIGNIFICANT  DATA  POINTS.  OPENS 
A  PILE  FOR  STORING  PLOTABLE  DATA.  AND  CALLS 
SUBROUTINE  PDPLOT.  WHICH  OBTAINS  A  PARAMETRIC 
FIT  TO  THE  DATA  AND  WRITES  DATA  AND  FIT  CURVES 
TO  A  FILE 

INTE0ER*4  NTRIAL.  NTC.  NTCX.NTCXO.  NTCMIN.  NTCMAX 


INTE0ER*4  NTRIAL.  NTC.  NTCX.NTCXO.  NTCMIN.  NTCMAX 
9IN8ERT  JIMB  JOD I >SUR ADS> I OLUS 

DIMENSION  NTC < 1 ) .  SBR(l).  PFAPCT ( 20 ), PDPCT (20).  S 
DATA  NTCMIN. I0UT/000010.  10/ 

C 

C  READ  PLOT  FILENAME  AND  OPEN  IT 
C 

WRITE  (LUOUT. 100) 

100  FORMAT  (////IX. 'PLOT  FILE  GENERATOR  ') 

READ  (LUIN.  llOHlFNAM(I).  1-1,3) 

110  FORMAT (3A2) 

READ (LUIN,  *) (PFAPCT ( J), J«l, NTHR) 

CALL  SEARCH  (3.  IFNAM,  IOUT. *800 > 

IWRIT-IOUT+4 

WRITE  (LUOUT,  120)  ( IFNAM( I ) .  1-1,  3) 

120  FORMAT  (//6X,  'PLOT  FILE  '.3A2,  '  OPENED' 

*//6X. 'FIT  PARAMETERS:  K1  K2  K3 

C 

C  LOOP  ON  THRESHOLDS 

C 

NTCMAX-NTR I AL-NTCMIN 
XTRIAL-100.  /FLOAT (NTRIAL) 

DO  200  1*1,  NTHR 
IJ-I-NTHR 
NPTS-O 
NTC X 0-0 
PFAT —PFAPCT ( I ) 

NTCMN-PFAT*.  01*2.  0*FL0AT( NTRIAL) 

NTCMN-MAXO< NTCMIN.  NTCMN) 

DO  130  J-l.NSBR 
I J— I J»NTHR 
NTCX-NTC(IJ) 

IF(.  NOT.  (  (NTCX.  OE.  NTCMN).  AND. 

1  (NTCX.  LE.  NTCMAX) ) )  GO  TO  130 

IF ( .  NOT.  (NTCX.  OE.  NTCXO)  >  GO  TO  130 

IF(.  NOT.  (NTCX.  NE.  NTCXO)  )  GO  TO  130 

NPTS-NPTS+1 

130  CONTINUE 

PDT-FLOAT ( NTC  X ) *XTR I AL 

PDPCT(NPTS)-100.  #(PDT— PFAT) / ( 100.  -PFAT) 
SDB(NPT8)-8BR(J) 

NTCXO-NTCX 
130  CONTINUE 

C 

C  FIT  DATA  AND  WRITE  TO  PLOT  FILE 

C 

CALL  PDPLOT  (PDPCT,  SDB. NPTS.  PFAT,  IWRIT) 

C 

200  CONTINUE 
C 

C  CLOSE  PLOT  FILE 
C 

CALL  SEARCH  (4.  IFNAM,  IOUT.  *900) 

WRITE  (LUOUT, 230)  ( IFNAM( I ),  I-l, 3) 

230  FORMAT  (//6X.  'PLOT  FILE  '.3A2,  '  CLOSED') 

C 

RETURN 

C 

900  STOP  800 

900  STOP  900 

END 


PDPCT(20).  SDB (20).  IFNAM(3) 


PFA(PCT) ' ) 
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A.  1.22  Generate  Plottable  Flies  (PDPLOT) 


Calling  Statement:  Call  PDPLOT  fPDPCT,  SDB,  NPTS,  PFAPCT,  IWRIT) 
Argument  List  Definitions; 

PDFCT  -  distorted  threshold  crossing  probabilities  in  percent 
SDB  “  corresponding  signal-tobackground  ratios  in  dB* 

NPTS  -  length  of  PDPCT  and  SDB 

PFAPCT  *  threshold  crossing  probability  for  no  signal 
IWRIT  -  Fortran  logical  unit  for  vriting  plottable  files. 

Other  Modules  Called:  QPDFIT ,  SBRVPD,  PDVSBR 

Description: 

PDPLOT  generates  plottable  files  of  PQ  vs  signal-tobackground 
data  and  parametric  fit  curves.  The  first  step  is  the  determination  of  the 
parameter  values  which  fit  the  distorted  PQ  data.  This  is  done  by  a  call  to 
QPDFIT.  Next,  these  parameter  vlaues  are  used  to  determine  the  signal-to¬ 
background  ratio  required  for  distorted  probabilities  of  1Z  and  99Z.  This 
is  done  by  a  call  to  SBRVPD.  A  smooth  curve  is  then  generated  by  a  call  to 
PDVSBR,  which  returns  smooth  fit  distorted  probabilities  at  each  of  25 
signal -to-background  ratios  spaced  evenly  in  dB  between  the  1Z  and  99Z 
values.  Finally,  the  distortion  is  removed  from  the  PQ  data  and  curve  fits, 
and  the  data  and  smooth  curve  pairs  written  to  a  file  for  off-line  plotting. 
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SUBROUTINE  PDPLOT  (PDPCT,  SDB> NPTS,  PFAPCT, I WRIT) 


C 

C 

C 

C 


GENERATE  PLOT ABLE  FILES  OF 

PD  VS.  SBR  DATA  AND  PARAMETRIC  FIT 


•INSERT  JIMB  JQD I >SUR ADS> I OLUS 

DIMENSION  PDPCT ( 1 ) ,  SDB  < 1 >  #  PDFITt 100) , SDBFIT ( 100) , 
«  P(3)» IFNAM<3)» PD199(2) 

DATA  IOUT,  XEOF, NFIT/ 10,  1.  E30.  23/ 

DATA  PD199/1.  , 99.  / 


C 

C 

C 


113 

C 

C 

C 

C 

C 

C 

C 

C 

C 


10 


20 

C 

C 

C 


IF(NPTS.  LT.  2)  00  TO  40 
FIT  THE  DATA 

CALL  GPDFIT<PDPCT#  SDB* NPTS#  P) 


IWRIT-IOUT+4 

WRITE  (LUOUT, 113)  (Pd), 

FORMAT </2 IX,  4010.  3) 


1*1,  3),  PFAPCT 


FIND  IX  AND  99X  SBR) 

CALL  SBRVPD(PD199, SDBFIT*  2#  P  > 

GENERATE  SMOOTH  FIT 

CALL  PDVSBR<PDFIT,  SDBFIT, -NFIT,  P) 

REMOVE  PFA  DISTORTION  FROM  DATA  AND  FIT 
DO  10  1*1, NPTS 

PDPCT ( I >*PDPCT ( I )•( ( 100.  -PFAPCT)/ 100.  >+PFAPCT 
CONTINUE 
DO  20  1*1, NFIT 

PDFIT ( I )*PDFIT < I )*<  < 100.  -PFAPCT)/100.  )+PFAPCT 
CONTINUE 


WRITE  DATA  AND  FIT 


120 


WRITEdWRIT,  120)((SDB(I),  PDPCT(I)  >, 
WRITEdWRIT,  120)  (XEOF.  XEOF) 

WRITEdWRIT,  120)  (  (SDBFIT  (I ) ,  PDFIT  (I  )  ) ,  I 
WRITEdWRIT.  120)  XEOF.  XEOF 
F0RMAT(2E12.  4) 


l.NPTS) 

1, NFIT) 


C 

C 

40  RETURN 
END 


A. 1.23 


in  Che 


where 

ratio 

[2]. 


Parametric  Curve  Fit  (QPDFIT) 

Calling  Statement:  Call  QPDFIT  (PD,  SDB,  N,  P) 

Argument  List  Definitions: 

PD  »  detection  probability  in  percent 
SDB  ■  signal-to-background  ratio  in  dB 
N  ■  length  of  PD  and  SDB 

P  ■  vector  of  3  coefficients  of  parametric  fit  determined 
by  QPDFIT. 

Description : 

QPDFIT  computes  approximate  values  for  parameters  P^,  P^,  and  P^ 
expression 

PD-100  x  f  1  -  [1  +  (P^)^] 


is  detection  probability  in  percent  and  S  ■  s ignal- to -background 
SDB 

log  !0  (^o“).  The  method  of  calculation  is  discussed  in  Reference 
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oo  ooo  ooo  ooo  non  ooo  oon  ooo  ooooooooooo 


SOB ROUTINE  QPDFXT  (PD# SDB#N#P> 

COMPUTES  APPROXIMATE  VALUES  FOR  PARAMETERS 
P1#P2#P3  IN  THE  EXPRESSION: 

PD-(1-(1+(P1*S)**P2)**-P3>*100 

WHERE:  PD-DETECTXON  PROBABILITY  (IN  PERCENT) 
S-S/N-AL0010(SDB/10.  0) 

SEE  TSC-W7-73 

DIMENSION  Y(2>#  Z(2) » W(200) # PD( 1 ) .  P (3) .  SDB(l) 
DIMENSION  PDD(200) 

REAL  MY#  MZ 
DO  100  I«1#N 

PDD(I)-PD(I>/100.  0 
CONTINUE 

COMPUTE  SLOPE  BE1WEEN  FIRST  TWO  POINTS 
DO  190  1-1#  2 

Y(I)-AL0010(-AL0010(1.  -PDD(I) >) 

CONTINUE 

MY-(Y(2)-Y(1))/(SDB(2)-SDB(1>> 

COMPUTE  SLOPE  BETWEEN  LAST  TWO  POINTS 

NM1— N— 1 
IZ-0 

DO  200  I-NM1,  N 
IZ-IZ+1 

Z(IZ>  —  AL0010(1.  -PDD(I)) 

CONTINUE 

HZ— (Z(2)-*Z( 1 )  ) /(SDB(N)-SDB (NM1 )  > 

COMPUTE  INITIAL  ESTIMATE  OF  K3 
P(3)tiZ/MY 

DISTORT  THE  PD  VALUES 
DO  229  I-l.N 

PDD(  I  >— 1.  -((1.  — PDD(  I )  )**(  1.  /P(3) ) ) 

CONTINUE 

COMPUTE  LARGE  PD  SLOPE  RATIO 

SRL1-AL0010(((1.  -PDD(N) )**(-!.  /P(3)>>-1.  ) 
SRL2-AL0010(((1.  -PDO(NMl) >**<-1.  /P(3>>)-1.  > 
8RL3-C-1.  /P(3)  >«*AI  0010(1.  -PDD(N)  > 

SRL4-( 1.  /P(3) )«AL0010( 1.  -PDD( NM1 >  > 

SRL- ( SRL1-SRL2 ) / ( SRL3+SRL4 ) 

COMPUTE  THE  SMALL  PD  SLOPE  RATIO 

SRS1«AL0G10( ((1.  -PDD(2) )**(-l.  /P(3>>>-1.  ) 
SRS2-AL0010(((1.-PDD(1> >**(-!.  /P(3)))-l.  ) 
SRB3-AL0©10( (-ALOOC 10.  0)/P(3) )*(AL0G10(1.  -PDD(2) ) ) ) 
SR84«AL0010(  (-ALOOUO.  0)  /P(3>  )»(AL0G10(  1.  -PDD(l)))) 
8R8- ( SRS 1 -SRS2 > / ( SRS3-SRS4 ) 

P(3>— P(3)*(SRL/8RS> 

INITIALIZE  COUNTERS  FOR  LINEAR  REGRESSION 

SUM-0 
SUMX-O.  0 
8UMY-0.  0 
SUM X 2-0.  0 
8UMXY-0. 0 

COMPUTE  W(X>  AND  PARAMETERS  FOR  A  STRAIGHT  LINE  FIT 


f 


ooo 


00  290  I-l.N 

W(I)«AL0010(((1.  01*PD< I ) )**<— 1.  /P<3) > )— 1.  > 

SUM-SUM* 1. 0 

SUMX«SUHX+SDB<I> 

8UMY«SUMY*W< I ) 

SUMX2-8UMX2+SDB ( I ) *80B ( I ) 

8UMXY-SUMXY+8DB ( I >  *W( I > 

CONTINUE 

D-3UH*SUMX2-SUHX*SUMX 

B- ( SUMX2*SU»1Y-SUMX*SUMXY ) /D 

A- ( SUMX Y*SUM-SUMX*SUHY  > /D 


COMPUTE  K1  AND  K2 


P<1)-10.  **<B/(10.  *A) ) 
P<2>-10.  *A 

RETURN 

END 


A.  1.24  SBR  vs.  PD  (SBRVPD) 

Calling  Statement:  Call  SBRVPD  (PDPCT,  SDB,  N,  P) 

Argument  List  Definition: 

PDPCT  “  detection  probability  in  percent 
SDB  ■  slgnal-to-baclcground  in  dB 
N  -  length  of  PDPCT  and  SDB 
P  ■  vector  of  coefficients 


where 


Description: 


Computes 


SDB  -  10  log. .  (S) 


■  t>[((i-v1/p0-1] 


PD  -  .01  x  PDPCT 


Pl’  P2’  P3  *  P(1)»  P(2)»  P(3) 


SDB  is  evaluated  at  each  of  N  values  of  PDPCT,  if  N  is  positive.  If  N  is 
negative,  SDB  is  evaluated  at  each  of  |n|  points  uniformly  spaced  between 
PDPCT  (1)  and  PDPCT  (2). 
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SUBROUTINE  SBRVPD < PDPCT. 8DB.  N. P ) 


COMPUTES: 


SOB* 10.  *AL0010(S(PD) ) 

WHERE* 

S-(i. /Pl>#<(<(l-PD>**(-l/P3>>-i>**(l/P2>> 

PD".  01*PDPCT 

PI.P2.  P3-P<1>.  P(2).  P(3) 


DIMENSION  P  <  3 }  #  POPCT ( 1 ) «  SDB ( 1 ) 

3(PX)-(1.  /P< 1 ))*<(( (1. 01*PX )**(-!.  /P<3) ) >—l.  )**(!.  /P(2) ) ) 


NI-IABS(N) 

IF ( .  MOT.  <N.  LT.  0)  >  00  TO  10 

PDSTEP-< PDPCT (2 > -PDPCT C 1 >  >/FLOAT<NI-l ) 

PD-PDPCT ( 1 ) — PDSTEP 

CONTINUE 

DO  20  1-1. NI 

IF(N.  OT.  0)  PD"PDPCT <  I  ) 

IF(N.  LT.  0)  PD-PD+PDSTEP 
PDPCT ( I )«PD 

SDK  I  >-10.  *AL0010(S(PD)  > 

CONTINUE 

RETURN 

END 


A.  1.25  PD  vs.  SBR  (PDVSBR) 


where 


Calling  Statement:  Call  PDVSBR  (PDPCT,  SDB,  N,  P) 
Argument  List  Definition: 

Same  as  for  SBRVPD 


Description: 

Computes 


PDPCT  -  100.  x 


C 


1 


S  -  io(SDB/10> 

Pl’  P2*  P3  *  P(1)*  P(2)*  P(3) 


PDPCT  is  evaluated  at  each  of  N  values  of  SDB,  if  N  is  positive.  If  N  is 
negative,  PDPCT  is  evaluated  at  each  of  |n|  points  uniformly  spaced  between 
SDB(l)  and  SDB(2). 


ouuugou 


SUBROUTINE  PDVSBR (PDPCT,  SOB« N,  P ) 


10 


30 


COMPUTES: 

PDPCT-100.  *<1.  -Cl.  +<P1*S)**P2>**<-P3> ) 
WHERE: 

8*10.  **8DB/10. 

PI,  P2»  P3»P<1>,P<2>.P<3> 

DZMEN8Z0N  P<3>, POPCT( 1 >.  S0B(1 > 

8P(X)«10.  **<X/10.  ) 

PD(S)-1.  -Cl.  ♦<P(1>*S>**P<2>  >**(-P<3>  > 
NI-IABS(N) 

IF ( .  NOT.  (N.  LT.  0))  00  TO  10 
8DB8TP*(SPB<2>— 8DB( 1 ) > /FLOAT (NI—l ) 
8DBL0C-SDB ( 1 >-SDB8TP 
CONTINUE 
DO  20  I«i«  NI 

IF(N.  OT.  0)  8DBL0C-3DB< I ) 

IF (N.  LT.  0)  SDBL0C«SD8L0C+8DB8TP 
8DB  <  I  > "8DBL0C 

PDPCT ( I ) «100.  »PD ( 8P ( SDBLOC ) > 

CONTINUE 

RETURN 

END 
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EXAMPLE  RUN 


A.  2 

An  example  run  of  the  simulation  program  la  included  bars  llluaerating 
tha  calculation  of  va.  SBR  In  distributed  rain  cluttar  (Raleigh  modulated 
log-normal  clutter  cross-section,  3  dB  log-s tender d  deviation),  for -a  alngle 
pulse. 

Figure  A.  2  Hats  all  Input  parameters  which  must  be  supplied  for  the 
simulation.  Figure  A. 3  lists  those  particular  values  used  in  the  example  run. 
Figure  A. 4  shows  the  output  file  CFBODT(LIJOOT) .  The  plot  file  is  not  shown. 
Table  A.l  tabulates  the  allowable  range  of  values  for  all  program  inputs. 


w 


INPUT  PARAMETER  HEAD  BY 


NTRIAL,  NPDI,  M,  IPLOT,  NTRLPS  .  DETDFX 

ITYPE,  NTHR,  (THR(I) ,  1-1,  NTHR) ,  NSBR,  (SBR(I),  1-1,  NSBR)  DETDFX 

If  ITYPE  -  1: 

ISEEDN  NOISE 

If  ITYPE  -  2: 

ISEEDL  LGN08M 

SDB  LGNORM 

If  ITYPE  -  3: 

ISEEDN  WEIBOL 

C  WEIBOL 

If  ITYPE  -  4s 

ITYPE  (local  variable  in  CORBKG)  CORBKG 

ISEEDN  NOISE 

If  ITYPE  -  2: 

ISEEDL  LGNORM 

SDB  LGNORM 


If  ITYPE  -  3: 

ISEEDW 

C 

ISEEDT 

ISWRL 

If  IPLOT  -  1: 

IFNAM(I),  I  -  1,  3 
PFAPCT(I) ,1-1,  NTHR 


Figure  A. 2.  General  List  of  Inputs  Required 


A-74 


I!!!  If 


NTRIAL,  NPDI ,  M,  IPLOT,  NTRLPS: 

200,  1,  127,  1,  1 

ITYPR  JRHR,(THR(I) ,  I  -  1,  NTHR)  : 

4,  6,  32.1,  33.3,  34.2,  35,  35.6,  36.2 
NSBR,  (SBR(I) ,  I  -  1,  USSR): 

20,  -3,  -2,  -1,  0,  1,  2,...,  16 
ITYPKlocml  Co  CORBKG): 

2 

ISEBDH: 

0 

ISEEDL: 

0 

SDB: 

3.0 

ISEEDT: 

0 

ISWRL: 

0 

IFNAM(I) ,  I  -  1,  3 
PDR1 

PFAPCT(I) ,  1-1,  NTHR: 

1.14,  0.27,  0.080,  0.022,  0.0073,  0.0021 
NTRIAL,  NPDI,  M,  IPLOT,  NTRLPS: 

0,  0,  0,  0,  0 


Figure  A.3.  Example  Run  Inputs 


f 


\ 

t 
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Figure  A. 4  OUTPUT  FILE  (CFROUT) 


01/10/80  0:  0:49  SYSTEM  (USER  *16) 


DETECTION  PERFORMANCE  OF  D1CKE-FXX  CFAR 


N8CAN8  -  200  NRBXN8  -  1 

BACKGROUND  TYPE  -  4  NPDX  - 

THRESHOLD  TEST  LEVELS  : 

32.  1  33.  3  34.  2 

39.  6  36.  2 


I  PLOT  - 
1  M  -  127 


39.  0 


1 


SXGNAL-TO— BACKGROUND  RATIOS  <DB> 


3.00 

-2.  00 

-1.00 

1.00 

2.00 

3.00 

9.00 

6.00 

7.00 

9.  00 

10.0 

11.  0 

13.0 

14.0 

19.0 

0.  000 
4.  00 
8.  00 
12.  0 
16.0 


MAXIMAL  LENGTH  PN  SEQUENCE 
CODE  LENGTH  -  127 

SHIFT  RE0X8TER  LENGTH  -  7 

FEEDBACK  CONNECTIONS  «  76 

WAVEFORM: 


!• 

l. 

l.  : 

1. 

1. 

1. 

—1. 

mm  X  • 

—1. 

—  1. 

-1. 

—1. 

1. 

—1. 

-1. 

1* 

-l. 

-l.  : 

1. 

-1. 

-1. 

-1. 

-1. 

1. 

-1. 

1. 

-1. 

-1. 

-1. 

1. 

l. 

1.  -! 

.  -1. 

1. 

-1. 

-1. 

-1. 

1. 

-1. 

1. 

1. 

-1. 

-1. 

1. 

1. 

-1.  3 

.  -1. 

1. 

-1. 

-1. 

1. 

1. 

1. 

1. 

1. 

-1. 

1. 

-1. 

-1. 

-1.  1 

1. 

1. 

-1. 

-1. 

-1. 

1. 

-1. 

-1. 

1. 

-1. 

-1. 

1. 

-1. 

1.  3 

.  -1. 

1. 

-1. 

1. 

1. 

-1. 

1. 

1. 

1. 

1. 

-1. 

1. 

a  • 

-1. 

-1.  -3 

.  1. 

1. 

-1. 

1. 

-1. 

-1. 

1. 

-1. 

1. 

1. 

1. 

-1. 

6  • 

1. 

1.  -3 

.  -1. 

1. 

1. 

•1. 

-1. 

1. 

mm  1 

1. 

-1. 

— j. 
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Figure  A. 4  (continued) 


CORRELATED  BACKGROUND  GENERATOR 

SPATIAL  MODULATION  TYPE  *  2 

GAUSSIAN  NOISE  GENERATOR 

PR N  GENERATOR  SEED  «  1739483429 

REPRESENTATION  CODE-  -1 

LOG-NORMAL  CLUTTER  GENERATOR 

PRN  GENERATOR  SEED  -  298809973 

LOO.  -STD.  -DEV.  -  3.  0  08 

REPRESENTATION  CODE  -  1 

TARGET  AMPLITUDE  GENERATOR 

PRN  GENERATOR  SEFD  -  924922231 

SWERLINO  CASE  0  TARGET  MODEL 
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Figure  A. 4  (continued) 


SUMMARY  OF  THRESHOLD  CROSS I NO  PROBABILITIES <PCT.  > 
NO.  OF  TRIALS  -  200 


THRESHOLD  CONST AN I 

S/B  RATIO (DB) 

32.  1 

33.  3 

34.  2 

39.  0 

-3.  0 

4.  00000 

1.  00000 

0.  00000 

0.  OOOOO 

>2.0 

4.  00000 

1.  00000 

0.  90000 

0.  OOOOO 

-1.0 

4.  00000 

1.  00000 

0.  90000 

0.  OOOOO 

-0.  0 

.  4.  00000 

2.  00000 

0.  90000 

0.  OOOOO 

1.  0 

9.  00000 

2.  00000 

0.  90000 

0.  OOOOO 

2.0 

9.  00000 

2.  90000 

0.  90000 

0.  OOOOO 

3.  0 

9.  90000 

3.  90000 

1.  90000 

0.  OOOOO 

4.  0 

8.00000 

3.  90000 

1.  90000 

0.  90000 

9.0 

11.  00000 

3.  90000 

2.  00000 

0.  90000 

6.  0 

12.  00000 

6.  90000 

2.  90000 

1.  OOOOO 

7.  0 

18.00000 

8.  90000 

4.  90000 

1.  90000 

8.  0 

27.  00000 

14.  90000 

6.  OOOOO 

3.  OOOOO 

9.0 

39.  00000 

23.  00000 

10.  90000 

6.  90000 

10.  0 

99.  90000 

34.  90000 

20.  OOOOO 

10.  90000 

11.  0 

69.  90000 

90.  90000 

31.  90000 

18.  OOOOO 

12.  0 

84.  00000 

69.  90000 

48.  OOOOO 

31.  90000 

13.  0 

93.  00000 

79.  90000 

69.  90000 

49.  OOOOO 

14.  0 

97.  00000 

92.  90000 

83.  OOOOO 

67.  90000 

19.  0 

99.  00000 

97.  90000 

99.  OOOOO 

89.  OOOOO 

16.  0 

99.  90000 

99.  90000 

98.  OOOOO 

96.  OOOOO 
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Figure  A. 4  (continued) 


THRESHOLD  CONSTANI 

S/B  RATXO(DB) 

39.  6' 

36.  2  ' 

-3.  0 

0.  00000 

0.  00000 

-2.  0 

0.  00000 

0.  00000 

-1.  0 

0.  00000 

0.  00000 

-0.  0 

0.  00000 

0.  00000 

1.  0 

0.  00000 

0.  00000 

2.0 

0.  00000 

0.  00000 

3.  0 

0.  00000 

0.  00000 

4.  0 

0.  00000 

0.  00000 

9.  0 

0.  00000 

0.  00000 

6.  0 

0.  00000 

0.  00000 

7.0 

0.  00000 

0.  00000 

8.  0 

1.  90000 

0.  00000 

9.0 

2.  90000 

0.  90000 

10.  0 

9.  90000 

2.  OOOOO 

11.  0 

11.  00000 

6.  OOOOO 

12.  0 

22.00000 

11.  90000 

13.0 

34.00000 

21.  OOOOO 

14.0 

94.  90000 

41.  OOOOO 

19.  0 

79.  90000 

61.  90000 

16.  0 

91.  90000 

S3.  OOOOO 

PLOT  FILE  GENERATOR 


PLOT  FILE  PDR1 

OPENED 

FIT  PARAMETERS: 

K1 

K2 

K3 

0.  497E-02 

t.  73 

136. 

0.  190E-01 

2.  10 

22.  3 

0.  274E-01 

2.  94 

9.  60 

0.  287E-01 

2.  66 

3.  26 

0.  294E-01 

3.  06 

2.  49 

0.  279E-01 

3.  33 

1.99 

PLOT  FILE  PDR1 

CLOSED 
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PFA(PCT) 

1.  14 
0.  270 
0.  800E-01 
0. 220E—01 
0.  730E-02 
0.  210E-02 


■  <** 

r : 


>  w*i  mrnmm 


Table  A.l.  Allowable  Ranges  of  Input  Variables 


Variable 

Type 

Allowable  Range  of  Values 

NTR1AL 

double  length  Integer 

1  <  NT RIAL  <  maximum  positive  double 
Integer  for  machine 

NPDI 

Integer 

NPDI  >  1 

M 

integer 

2  <  M  <  256 

( (M-l) *2  +  1)NPDI  <  2 04 A 

IPLOT 

Integer 

any  single  precision  Integer  value 

NTRLPS 

Integer 

1  <  NTRLPS  <  NTRIAL 

ITYPE 

integer 

1,  2,  3,  4 

NTHR 

integer 

1  <  NTHR  <  20 

THR  (20) 

real 

any  single  precision  real  value 

NSBR 

integer 

1  <  NSBR  <  20 

SBR  (20) 

real 

any  single  precision  real  value 

ISEEDN 

double  length  integer 

any  non-negative  double  Integer 

ISEEDL 

double  length  Integer 

any  non-negative  double  Integer 

ISEEOW 

double  length  Integer 

any  non-negative  double  integer 

ISEEDT 

double  length  Integer 

any  non-negative  double  Integer 

SDB 

real 

any  single  precision  real  value 

C 

real 

any  positive  single  precision  real  value 

ISWRL 

Integer 

0,  1,  2,  3,  4 

IFNAM  (3) 

integer 

any  six-character  file  name  (2  char. /word) 

PFAPCT  (20] 

real 

0.0  <  PFAPCT  <  100.0 
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